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We  consider  an  M/M/2  queueing  system  with  removable  service 
stations  operating  under  steady- state  over  an  infinite  time  horizon.  We 
assume  that  the  number  of  operating  service  stations  can  be  adjusted  at 
customers'  arrival  and  service  completion  epochs  depending  on  the  number 
of  customers  in  the  system. 

First,  the  steady-state  behavior  of  the  system  is  studied  as  a 
birth  and  death  process.  The  joint  probability  mass  function  of  the 
number  of  customers  in  the  system  and  the  number  of  operating  service 
stations  is  obtained  as  well  as  other  system  characteristics  which  might 
be  incorporated  in  a total  expected  cost  function  to  optimize  the 
operation  of  the  system. 

Next,  the  busy  period  is  analyzed  using  the  properties  of  the 
classical  gambler's  ruin  problem.  Fundamental  information  of  system 
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behavior  such  as  the  probability  mass  function  of  the  number  of 
activations  or  removals  of  service  stations  during  the  busy  period  and 
the  Laplace  transform  of  the  probability  density  function  of  the  busy 
period  are  derived.  As  a special  case,  a number  of  existing  results 
such  as  the  Laplace  transforms  of  the  probability  density  functions  of 
the  busy  periods  for  the  ordinary  M/M/2  system,  the  M/M/1  system 
operation  under  the  N policy  and  the  ordinary  M/M/1  system,  are 
recovered.  Alternative  methodologies  are  also  provided  to  find  certain 
system  characteristics  of  particular  use  in  a cost  model. 

Finally,  a total  expected  cost  function  per  unit  time  is 
constructed  to  obtain  the  optimal  customer  levels  dictated  by  the  three 
decision  variables  at  which  adjustments  of  service  stations  are  made. 

A sensitivity  analysis  is  also  performed  for  changes  of  the  optimal 
customer  levels  along  with  changes  in  specific  cost  elements. 
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CHAPTER  1 


INTRODUCTION 


1 • 1 General  Description  of  Controllable  Queueing  Models 

Almost  all  the  fields  of  modern  society  have  faced  the  problem  of 
providing  huge  demands  for  simultaneous  service.  This  reality  makes 
waiting  lines  an  unavoidable  component  of  modern  life  style  leading  to 
accelerated  research  activities  in  waiting  line  systems.  The  general 
trend  of  research  has  been  to  create  mathematical  models  appropriately 
based  on  real  waiting  line  situations  and  to  derive  important  system 
characteristics  using  these  models. 

In  the  early  stages  of  queueing  theory,  the  basic  assumptions  for 
the  creation  of  mathematical  models  involved  fixed  system  parameters. 

In  other  words,  each  system  parameter  could  not  be  allowed  to  change 
over  time  or  system  state.  The  basic  assumptions  were  critical  factors 
not  only  to  describe  in  a simple  fashion  the  real  situations,  but  to  get 
desired  results  easily.  However,  although  the  results  obtained  could 
provide  information  to  predict  real  system  behavior,  there  seems  to  be  a 
large  gap  between  assumed  models  and  real  situations  because  the  basic 
assumptions  might  be  far  from  realistic.  This  motivates  a change  of 
research  direction  by  considering  models  which  could  afford  relaxations 
of  those  assumptions. 

The  controllable  queueing  models  which  allow  the  system  parameters 
to  be  changed  over  time  or  system  state  have  been  a major  part  of 
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research  in  queueing  theory  over  the  last  two  decades.  Numerous 
different  models  which  could  be  considered  as  control  models  have  been 
introduced  in  the  literature.  There  could  be  a nirmber  of  different 
criteria  to  classify  the  introduced  models.  However,  Crabil,  Gross  and 
Magazine[9]  and  J.Teghem  Jr. [25]  distinguish  four  categories  for  these 
models  depending  on  which  system  parameter  can  be  changed  over  time. 
These  categories  are  'control  of  the  number  of  servers  or  service 
stations',  'control  of  the  service  rate',  'control  of  the  admission  of 
the  arrival  customers'  and  'control  of  the  waiting  line  discipline'. 

Here  we  are  concerned  with  the  models  related  to  the  'control  of 
the  number  of  service  stations';  the  so  called  vacation  models.  The 
main  purpose  of  the  vacation  models  is  to  try  to  utilize  the  server's 
time.  There  could  be  partial  or  complete  removals  of  service  stations 
based  on  a predetermined  operating  policy.  Thus,  the  criteria  for 
activations  and  removals  of  the  service  stations  and  the  epochs  at  which 
those  adjustments  should  occur  are  the  most  important  assumptions  of  the 
operating  policy.  Review  points  or  decision  points  are  defined  on  the 
time  axis  as  those  epochs  at  which  the  decisions  are  made  on  whether  the 
adjustments  of  the  service  stations  should  occur  or  not.  According  to 
the  definition  of  review  points,  the  vacation  models  can  be  classified 
into  three  categories.  The  model  is  called  the  periodic  review  model  if 
the  time  intervals  between  any  two  successive  review  points  are  the  same 
and  determinations  of  the  review  points  are  independent  of  the  system 
states.  If  the  review  points  are  stochastically  determined  depending  on 
the  system  states,  then  the  model  is  called  the  stochastic  review  model. 
In  this  model,  the  review  points  are  not  necessarily  equally  spaced  in 
time.  The  last  type  is  the  continuous  review  model  in  which  the 
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adjustments  can  be  made  continuously  over  time.  The  customers'  arrival 
and  service  completion  epochs  are  one  of  the  possible  choices  as  review 
points  in  the  area  of  vacation  models.  The  criteria  for  adjustments  of 
service  stations  are  closely  related  to  the  decision  variables.  Based 
on  the  activation  criteria  of  removed  service  station,  three  different 
types  of  operating  policy  have  been  introduced.  The  first  one  is  called 
the  N policy  in  which  the  decision  variables  for  activations  of  service 
stations  are  the  number  of  customers  in  the  system  (see  Yadin  and 
Naor[27]).  Thus,  the  removed  service  station  can  only  be  activated  at  a 
time  when  the  number  of  customers  in  the  system  reaches  a specific 
level.  If  the  decision  variables  are  workloads  which  are  the  total 
amounts  of  existing  customers'  work  in  the  system  or  server's  backlog, 

D,  then  the  operating  policy  is  called  the  D policy  (see  Balachandran 
and  Tijms[l]  and  Boxma[6]).  Finally,  an  operating  policy  is  called  the 
T policy  if  service  stations  only  activate  T time  units  after  the  end  of 
the  last  busy  period  (see  Heyman[12]).  Furthermore,  according  to  the 
removal  criteria  of  activating  service  stations,  one  can  have  the 
exhaustive  policy  or  the  hysteretic  policy.  If  the  activated  service 
stations  can  only  be  removed  at  a time  when  the  service  stations  become 
idle,  then  one  has  the  exhaustive  policy.  The  other  case  is  the 
hysteretic  policy  in  which  service  stations  are  removed  even  when  they 
do  not  become  idle . 

The  general  approach  used  in  the  controllable  queueing  problems  is 
to  decide  on  the  optimal  operating  policy.  This  can  be  done  based  on 
minimizing  the  total  expected  cost  involved  in  system  operation. 

Usually,  the  cost  elements  consist  of  holding  cost,  operating  cost,  and 
switching  cost.  The  holding  cost  depends  on  the  number  of  customers  in 
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the  system.  The  operating  cost  is  related  to  the  number  of  operating 
service  stations.  The  adjustments  of  service  stations  contribute  to  the 
switching  cost. 


1 . 2 Literature  Review 

The  paper  of  Moder  and  Phillips  Jr. [19]  is  a typical  example  which 
shows  the  general  research  pattern  dealing  with  the  control  of  the 
number  of  seirvice  stations  in  the  beginning  years  of  controllable 
queueing  models.  During  this  time,  the  prime  interest  of  researchers 
was  to  derive  the  only  system  characteristics  which  could  be  related  to 
cost  elements  under  the  assumed  operating  rules.  Even  though  the 
reasonable  operating  rules  were  well  developed,  optimality  of  the 
assumed  operating  rules  was  not  discussed. 

More  recently  the  interests  of  researchers  in  the  area  of  removable 
servers  or  service  stations  have  changed.  The  developments  of  realistic 
operating  policies  such  as  the  N policy,  the  D policy  and  the  T policy, 
and  the  determinations  of  the  optimal  decision  variables  have  been  of 
major  interest.  However,  we  concentrate  on  those  papers  dealing  with 
the  N policy.  At  this  stage,  it  should  be  pointed  out  that  most  of  the 
work  in  the  area  of  N policy  are  concerned  with  models  with  a single 
removable  service  station.  Only  a few  papers  devote  themselves  to  the 
multi -removable  service  stations  operating  under  the  concept  of  the  N 
policy. 

The  next  section  is  divided  into  two  parts  according  to  the  number 
of  removable  service  stations  studied  in  the  literature.  The  first  part 
reviews  those  papers  which  treat  the  single  removable  service  station 
operating  under  an  N policy.  The  second  part  deals  with  those  papers 
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treating  the  multi -removable  service  station  models. 

1.2.1  Single  Removable  Service  Station  Operating  under  an  N policy 

The  concept  of  an  N policy  in  queueing  theory  was  first  introduced 
by  Yadin  and  Naor[27] . They  consider  an  M/G/1  system  in  which  the 
server  is  removed  at  the  end  of  a busy  period.  Then,  the  server  is 
activated  when  the  number  of  customers  in  the  system  reaches  R for  the 
first  time.  They  derive  the  expected  number  of  customers  in  the  system 
for  the  removable  M/G/1  system.  They  also  derive  the  optimal  decision 
variable  R which  minimizes  the  total  expected  cost  per  unit  time 
including  the  holding  costs,  the  switching  costs  and  certain  negative 
costs  (profits)  which  are  the  results  of  utilization  of  server's  idle 
time . 

Heyman[ll]  also  considers  an  M/G/1  queueing  system  with  removable 
service  stations.  His  cost  elements  consist  of  the  operating  cost,  the 
holding  cost  and  the  switching  cost.  He  investigates  the  optimal 
behaviors  of  the  system  under  the  average  cost  and  the  discounted  cost 
criterion  over  an  infinite  time  horizon.  For  the  average  cost 
criterion,  he  derives  explicitly  the  optimal  customer  level  at  which  the 
service  station  should  be  opened.  For  the  discounted  cost  criterion,  he 
proves  the  same  general  properties  of  the  optimal  operating  N-policy  as 
those  of  the  average  cost  criterion.  However,  because  of  the 
difficulties  in  finding  the  optimal  customer  level  explicitly,  he 
provides  an  algorithm  to  determine  it  indirectly. 

Sobel[22]  studies  a GI/G/1  removable  service  station  with  an  N 
policy  using  the  average  cost  criterion  over  an  infinite  time  horizon. 

He  shows  that  Heyman's  results  on  the  form  of  the  optimal  operating 
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policy  are  valid  for  his  model. 

Kimura[14]  applies  diffusion  approximation  method  to  Heyman[ll]'s 
discounted  cost  model.  He  uses  diffusion  processes  to  approximate  the 
number  of  customers  in  the  system  by  a continuous  variable.  He  shows 
numerically  that  his  approach  can  generate  the  Heyman's  results  under 
heavy  traffic  condition. 

Several  interesting  extensions  of  the  above  models  can  also  be 
found  in  the  literature.  But  we  only  introduce  here  a brief  description 
of  these  models.  A more  detailed  review  appears  in  the  paper  of  Teghem, 
Jr. [25]. 

Some  researchers  have  studied  the  single  removable  server  queueing 
system  operating  under  an  N policy  for  a finite  system  capacity.  For 
example,  one  can  cite  Hersh  and  Brosh  (see  Teghem,  Jr. [25])  for  the 
M/M/1  system,  and  Lolis  and  Teghem  (see  Teghem,  Jr. [25])  for  the  M/G/1 
system.  In  these  types  of  models,  the  holding  costs  are  usually 
replaced  by  the  penalty  costs  incurred  for  lost  customers. 

Others  have  extended  the  basic  models  using  priority  queueing 
system.  Bell [2, 3]  considers  an  M/G/1  priority  queue  with  removable 
server  using  an  infinite  horizon  average  cost  as  the  criterion.  He 
assumes  several  priority  classes  of  customers,  each  class  of  customer 
having  the  different  holding  cost  and  the  same  service  time 
distribution.  Tijms  (see  Teghem,  Jr. [25])  studies  an  M/G/1  non- 
preemptive  priority  queueing  system  with  a removable  server  and  two 
priority  classes  of  customers.  He  assumes  that  each  class  of  customers 
has  a different  holding  cost  and  a different  service  time  distribution. 
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1.2.2  Multi -removable  Service  Stations 

Magazine [16]  considers  an  M/M/s  queueing  system  with  removable 
service  stations  under  periodic  review.  The  number  of  customers  in  the 
system  is  observed  at  equally  spaced  review  points  in  time.  His 
decision  variables  are  the  number  of  servers.  Decisions  may  be  made  at 
each  review  point.  If  the  observed  number  of  customers  in  the  system  at 
each  review  point  falls  between  two  predetermined  values,  then  a 
decision  is  made  on  how  many  servers  should  be  activated  or  removed.  If 
it  falls  outside  the  range  of  two  values,  no  action  is  taken.  Using 
dynamic  programming,  he  shows  the  existence  of  a control  limit  policy 
for  the  average  cost  criterion  over  a finite  and  an  infinite  time 
horizon.  In  other  words,  there  exists  a unique  number  of  customers  in 
the  system,  say  j*,  such  that  in  order  to  activate  the  service  station 
the  observed  number  of  customers  in  the  system  should  be  greater  than 
j*.  He  also  proves  the  control  limit  policy  is  optimal  for  the  infinite 
horizon  model  and  the  average  cost  criterion. 

Magazine [17]  extends  his  previous  model  under  assumption  of  a 
firiite  system  capacity.  He  also  derives  the  same  conclusions  as  those 
of  his  previous  model. 

Huang,  Brumelle,  Sawaki  and  Vertinsky[13]  consider  a similar  model 
to  the  one  outlined  by  Magazine [ 17 ] , with  some  modifications.  They 
assume  infinite  system  capacity  and  generalized  holding  costs  under  a 
discounted  cost  criterion  over  a finite  time  horizon.  They  permit 
different  holding  costs  which  are  functions  of  the  number  of  customers 
in  the  system  at  the  end  of  each  period.  They  also  assume  that  the 
holding  cost  functions  are  convex.  They  conclude  that  the  optimal 
operating  policy  has  the  form  of  a control  limit  policy. 


8 


Szarkowicz  and  Knowles [23]  investigate  the  model  studied  by  Huang 
et  al.[13]  under  less  restrictive  assumptions.  But  they  derive  the 
same  conclusions  as  those  of  Huang  et  al. 

In  the  paper  by  McGill [18],  an  M/M/s  and  a GI/M/s  queueing  system 
are  investigated  under  stochastic  review.  He  assumes  that  the  time 
intervals  between  any  two  successive  review  points  are  random  variables 
which  depend  on  the  system  states.  Like  Magazine [ 17 ] , his  decision 
variables  are  also  the  number  of  servers  to  be  activated  or  removed 
based  on  the  number  of  customers  in  the  system.  He  develops  a framework 
for  determination  of  the  optimal  decision  variables  under  total  expected 
discounted  and  undiscounted  costs  criteria  for  the  M/M/s  system.  He 
derives  the  explicit  optimal  control  policy  for  the  case  of  no  switching 
cost . 

Bell [4]  treats  an  M/M/s  model  studied  by  McGill.  The  number  of 
servers  working  can  be  adjusted  at  customers'  arrival  or  service 
completion  epochs.  He  defines  an  efficient  policy  as  an  operating 
policy  which  never  allows  more  active  service  stations  than  the  number 
of  customers  presented  in  the  system,  otherwise  an  operating  policy  is 
inefficient.  He  shows  numerically  that  the  optimal  operating  policy  may 
not  necessarily  be  an  efficient  policy. 

Winston [26]  examines  a removable  M/M/s  system  in  which  the 
customers'  arrival  rates  depend  on  the  current  system  state.  His  review 
points  are  also  the  customers'  arrival  and  service  completion  epochs, 
and  his  decision  variables  are  the  number  of  servers.  Unlike  other 
models , he  assumes  that  the  removal  or  the  activation  of  service 
stations  incurs  no  costs.  But  he  generalizes  the  operating  costs  and 
the  holding  costs.  His  holding  costs  are  a function  of  the  number  of 


9 


customers  in  the  system.  His  operating  costs  are  also  a function  of  the 
number  of  active  service  stations.  He  derives  optimal  conditions  that 
ensure  that  the  number  of  servers  in  operation  is  a non- decreasing 
function  of  the  number  of  customers  in  the  system. 

More  recently  Bell [5]  tries  to  apply  the  concept  of  an  N policy  to 
the  multi -removable  service  station  case.  He  considers  an  M/M/2 
queueing  system  where  the  customers'  arrival  and  service  completion 
epochs  are  the  review  points.  He  defines  the  (Sq  , ,R2)  policy 

denoting  the  number  of  customers  in  the  system  when  the  number  of 
opening  servers  should  be  adjusted  downward  to  0,1  and  upward  to  1,2, 
respectively.  He  allows  the  possibility  that  the  number  of  opening 
service  stations  is  never. adjusted  downward  to  0 or  1,  denoting  by  Sq  - 
-1,  or  - -1.  Without  deriving  the  system  characteristics  explicitly, 
he  shows  that  the  optimal  operating  policy  is  an  efficient  policy  if  the 
operating  costs  are  sufficiently  high.  If  the  operating  costs  are 
sufficiently  close  to  zero  compared  to  the  holding  costs  and  switching 
costs,  an  inefficient  policy  is  optimal.  He  also  proves  general 
properties  of  the  optimal  operating  policy.  For  the  efficient  policy, 
the  optimal  policy  has  the  form  Sq  = 0,  1 < < R2  For  the 

inefficient  policy,  the  optimal  operating  policy  may  be  one  of  following 
cases:  (1)  leave  both  servers  opened  at  all  times,  denoted  by  Sq  - - 

-1,  (2)  leave  at  least  one  server  opened  at  all  times,  denoted  by  Sq  - 
•1.  > 0,  Sj^  < R2,  or  (3)  always  close  both  servers,  denoted  by  Sq  - 

S 2^  “ 0 , R2  ^ R^^ . 

Even  though  most  researchers  have  concentrated  on  deriving  the 
optimal  operating  policy  rather  than  analyzing  the  system  behavior  in 
the  area  of  multi -removable  service  stations,  no  one  has  attempted  to 
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derive  it  explicitly  because  of  the  complexity  of  the  problem. 

1 . 3 Problem  Statement 

In  this  research,  we  consider  an  M/M/2  queueing  system  with 
removable  service  stations  under  steady- state  over  an  infinite  time 
horizon.  The  review  points  of  the  system  are  the  customers'  arrival  and 
service  completion  epochs.  The  decisions  on  the  adjustment  of  the 
service  stations  are  made  based  on  the  number  of  customers  in  the 
system.  Thus,  the  operating  policy  of  our  model  should  be  considered  a 
type  of  the  N policy.  The  decision  variables  of  our  model  are  several 
levels  of  customers  in  the  system  at  which  the  adjustments  of  the 
service  stations  should  occur. 

1.3.1  Operating  Policy  of  the  System 

Consider  an  M/M/2  steady-state  queueing  process.  The  customers' 
interarrival  times  are  identically,  independently  and  exponentially 
distributed.  The  customers'  service  times  at  both  service  stations  are 
also  identically,  independently  and  exponentially  distributed.  The 
customers'  interarrival  times  and  the  customers'  service  times  are 
independently  distributed.  The  special  feature  of  the  operating  policy 
of  this  queueing  system  is  that  the  number  of  operating  service  stations 
can  always  be  adjusted  based  on  the  number  of  customers  present  in  the 
system.  Thus,  the  number  of  customers  in  the  system  is  monitored  at 
every  new  customer's  arrival  and  service  completion  epochs.  Whenever 
there  are  no  customers  in  the  system,  both  service  stations  remain 
inoperative  temporarily  and  may  not  be  reactivated  until  certain 
conditions  are  satisfied.  Suppose  both  service  stations  are  removed 
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from  the  system.  When  the  number  of  customers  waiting  for  service 
reaches  a specified  level,  denoted  by  N which  is  a positive  integer,  one 
of  the  two  service  stations  removed  from  the  system  is  activated 
immediately.  At  a later  time  when  the  number  of  customers  in  the  system 
increases  to  another  specified  level,  say  M where  0 < N < M,  then  the 
remaining  service  station  removed  from  the  system  is  also  activated 
immediately.  However,  if  the  number  of  customers  in  the  system 
decreases  to  K where  1 < K < N while  both  service  stations  are 
operating,  the  service  station  just  finishing  a service  will  be  removed 
at  that  instant.  Furthermore,  if  the  number  of  customers  in  the  system 
reaches  zero  while  one  service  station  is  operating,  that  service 
station  should  also  be  removed  until  the  above  reactivating  conditions 
are  satisfied. 

We  denote  the  operating  policy  of  our  model  as  the  (0,K,N,M) 
policy.  This  operating  policy  is  an  efficient  policy  according  to 
Bell[4]'s  definition  because  the  number  of  active  service  stations  is 
always  less  than  the  number  of  customers  present  in  the  system. 
Furthermore,  no  service  stations  are  allowed  to  operate  when  the  system 
becomes  empty.  Thus,  this  operating  policy  can  be  considered  as  a type 
of  the  exhaustive  policy. 

1.3.2  Special  Cases  of  the  Model 

We  assumed  that  the  niomber  of  operating  service  station  could  be 
adjusted  based  on  the  number  of  customers  in  the  system.  Hence,  unknown 
values  of  the  decision  variables  K,  N and  M are  crucial  criteria  to 
determine  the  system  operation.  This  implies  that  each  value  of  the 
decision  variables  specifies  a different  operating  policy  of  the  system. 


12 


Now,  we  consider  two  typical  special  cases  of  the  (0,K,N,M)  operating 
policy  by  assigning  specific  values  to  the  decision  variables. 

g.ase  1 : Suppose  that  K=l,  N=1  and  M=2 , denoted  by  (0,1, 1,2)  policy. 

Then,  whenever  there  is  a customer  in  the  system,  the  number  of 
operating  service  stations  is  one.  If  the  number  of  customers  in  the 
system  is  greater  than  or  equal  to  two,  both  service  stations  are 
operating  simultaneously.  When  the  system  becomes  empty,  no  service 
stations  are  operating.  Hence,  it  is  easy  to  see  that  the  M/M/2 
queueing  system  operating  under  the  (0,1, 1,2)  policy  is  equivalent  to 
the  ordinary  M/M/2  queueing  system. 

2 : Suppose  that  M-«,  denoted  by  (0,K,N,=o)  policy.  One  of  the  two 

service  stations  removed  from  the  system  starts  operating  when  the 
number  of  customers  in  the  system  reaches  N.  Since  M=«>,  the  remaining 
service  station  removed  from  the  system  cannot  start  operating  at  all. 
Hence,  the  operating  service  station  can  only  be  removed  when  there  are 
no  customers  in  the  system.  In  this  case,  the  behavior  of  the  system  is 
independent  of  the  value  of  K.  Thus,  the  behavior  of  the  M/M/2  queueing 
system  operating  under  the  (0,K,N,«)  policy  becomes  that  of  the  M/M/1 
queueing  system  operating  under  the  N policy  which  is  studied  by  Yadin 
and  Naor [27] . 

We  thus  conclude  that  the  M/M/2  queueing  system  operating  under  the 
(0,K,N,M)  policy  is  considered  as  a generalized  model  of  the  ordinary 
M/M/2  queueing  system  and  the  M/M/1  queueing  system  operating  under  the 
N policy. 
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1.3.3  Applicable  Examples  of  the  Model 

Our  model  could  be  applicable  to  approximate  a number  of  real 
situations  such  as  banks,  post  offices,  production  planning, 
manufacturing,  communication  and  so  on.  We  elaborate  on  two  examples. 

First,  we  give  an  example  in  manufacturing  area. 

Example  1 : Consider  a manufacturing  plant  where  its  major  tasks  are 

assembly,  inspection  and  rework.  Suppose  that  a large  number  of  workers 
assemble  the  same  type  of  expensive  items  individually  at  the  assembly 
department.  All  assembled  items  are  tested  at  the  inspection 
department.  Whenever  inspectors  find  a defective  item,  they  send  it  to 
the  rework  department  one  by  one.  Interarrival  times  of  defective  items 
at  the  rework  department  are  exponentially  distributed.  Each  defective 
item  could  be  corrected  by  additional  work  or  could  be  disassembled  to 
reutilize  some  components  of  defective  items  at  the  rework  department. 
However,  all  the  defective  items  arriving  at  the  rework  department  are 
kept  in  special  type  of  containers  in  order  to  protect  them  from  further 
damage.  The  corrections  and  the  disassembly  for  defective  items  can  be 
done  effectively  by  ones  who  are  familiar  with  the  assembling  process . 
Approximately,  the  correction  or  the  disassembly  time  for  an  defective 
item  is  exponentially  distributed.  Workloads  at  the  rework  department 
vary  depending  on  the  number  of  defective  items  arrived.  It  is  not  a 
good  strategy  to  keep  the  workers  always  staying  at  the  rework 
department,  because  there  could  be  no  defective  items  to  be  reworked. 
Thus,  to  utilize  the  reworking  workers'  times,  the  rework  department 
should  be  inoperative  under  certain  conditions.  This  circumstance 
requires  that  an  operating  rule  of  the  rework  department  has  to  be 
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Figure  1.1  The  flow  diagram  of  assembling  process 
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related  to  the  number  of  defective  items  arrived.  In  order  to  analyze 
this  situation,  without  loss  of  generality,  assume  that  two  of  the 
assembly  workers  are  designated  to  be  part  time  crew  at  the  rework 
department.  Their  tasks  are  changed  alternately  depending  on  where  they 
work,  the  rework  or  the  assembly.  Whenever  there  are  no  defective  items 
arriving  at  the  rework  department,  the  rework  department  is  inoperative 
temporarily.  Thus,  two  designated  workers  perform  the  assembling  tasks 
like  other  assembly  workers  at  the  assembly  department.  If  the  number 
of  defective  items  arriving  at  the  rework  department  reaches  N,  then  one 
of  the  two  designated  workers  stops  the  assembling  task  and  immediately 
starts  the  reworking  task  at  the  rework  department.  After  that,  if  the 
number  of  arriving  defective  items  increases  to  M which  is  greater  than 
N,  the  other  designated  worker  also  joins  the  rework  department. 

However,  if  the  number  of  defective  items  decreases  to  K where  K < N 
while  both  workers  are  at  the  rework  department,  one  of  the  workers 
stops  the  reworking  task  and  returns  to  the  assembly  department  in  order 
to  perform  the  assembling  task.  The  remaining  worker  closes  the  rework 
department  and  returns  to  the  assembly  department  in  order  to  perform 
the  assembling  task  if  there  are  no  defective  items  at  the  rework 
department.  Otherwise,  the  remaining  worker  performs  reworking  task  for 
the  defective  items.  The  removals  and  the  additions  of  two  designated 
workers  at  the  assembly  department  do  not  significantly  change  the 
defective  rate  of  assembling  items.  This  operating  policy  of  the  rework 
department  can  not  only  utilize  the  worker's  time,  but  also  assign 
balancing  workloads  to  each  worker.  However,  this  operating  policy 
requires  switching  costs  of  the  reworking  machines  and  the  holding  costs 
of  defective  items  to  keep  in  the  containers  at  the  rework  department. 
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The  optimal  operating  policy  at  the  rework  department  has  to  be  selected 
based  on  minimizing  the  total  expected  costs. 

Another  example  is  related  to  a communication  problem. 

Example  2 : Suppose  that  the  main  office  of  a sales  company  receives 

important  marketing  information  or  messages  from  several  sources. 
Arrivals  of  the  messages  at  the  main  office  are  assumed  to  follow  a 
Poisson  stream.  All  messages  received  are  to  be  sent  to  a branch  office 
in  a foreign  country  by  a commercial  satellite  system.  The  company  can 
use  two  access  lines  via  the  satellite  to  send  the  messages.  The  access 
lines  which  can  be  considered  as  the  service  stations  are  available 
whenever  the  company  wants  to  send  the  messages.  Transmission  time  of 
each  message  via  the  satellite  is  exponentially  distributed.  But 
whenever  the  company  uses  any  available  access  lines,  the  company  has  to 
pay  more  charges  for  a first  service  time  unit  and  less  for  subsequent 
service  time  units.  Thus,  a turn-on  cost  is  incurred  in  the  use  of  the 
access  line  which  is  the  excess  charge  for  the  first  service  time  unit. 
On  the  one  hand,  it  may  be  economical  to  turn-on  both  access  lines 
continuously  and  to  send  the  messages  whenever  the  main  office  receives 
them.  Thus,  the  company  can  avoid  the  turn-on  costs  of  access  lines 
even  though  unnecessary  operating  costs  of  the  access  line  are  incurred 
when  there  are  no  messages  to  be  sent.  On  the  other  hand,  the  company 
may  experience  a certain  type  of  loss  in  using  the  access  lines 
selectively  because  the  messages  received  are  kept  for  a longer  time 
without  sending  them  to  the  branch  office.  Therefore,  contribution  of 
the  turn-on  costs,  the  operating  costs  of  the  access  lines  and  the 
holding  costs  of  the  messages  received  would  determine  the  optimal 
operating  policy  of  access  lines.  The  optimal  decision  variables  which 
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ths  total  6xp8ct6d.  cost  involvsd  at6  th©  levels  of  numbet  of 
messages  received  at  which  access  lines  should  be  turned-on  or  turned- 
off.  Therefore,  it  is  reasonable  to  opt  for  the  following  operating 
policy  for  the  access  lines:  The  main  office  keeps  the  messages 

received  without  sending  them  to  the  branch  office  until  their  number 
reaches  a certain  level,  denoted  by  N.  If  that  number  reaches  N,  one  of 
the  available  access  lines  is  used  and  messages  are  sent  in  the  order  of 
arrivals  until  they  are  all  sent  to  the  branch  office.  If  that  number 
reaches  a level  M > N,  the  other  access  line  starts  operating  until  the 
number  of  remaining  messages  reaches  a level  K < N at  which  time  one  of 
the  operating  access  lines  is  turned  off.  One  is  then  interested  in  the 
determination  of  the  optimal  decision  variables  (K*,N*,M*)  to  minimize 
the  total  expected  cost  per  unit  time. 

1-3.4  Basic  Assumptions  of  System  Operation 

To  formulate  mathematically  the  problem,  the  following  assumptions 
are  made . 

(1)  Customers'  arrivals  follow  a Poisson  process  with  intensity  A. 

In  other  words,  the  random  variables  where  i-1,2,3,---  which 
represent  the  interarrival  times  between  the  (i-l)th  and  the  ith 
customers  joining  the  system  are  identically,  independently  and 
exponentially  distributed  with  E[A^^]  = 1/A  for  all  i.  The  common 
probability  density  function  of  A^  is  denoted  by 

f^(t)  = Ae'^^  A > 0,  t > 0. 

(2)  The  random  variables  S^'s  where  i=l,2,3, • • • which  represent  the 
ith  customer's  service  time  are  identically,  independently  and 
exponentially  distributed  with  E[S^]  = l//i  for  all  i.  The  common 
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probability  density  function  of  is  denoted  by 
fgCt)  = M ^ 0,  t > 0. 

(3)  The  ' s and  the  S^'s  are  independently  distributed. 

(4)  Both  service  stations  are  indistinguishable. 

(5)  Arriving  customers  form  a single  waiting  line. 

(6)  The  service  stations  are  not  subject  to  break-downs. 

(7)  Customers  are  served  in  order  of  their  arrivals 
(First-come,  First-served). 

(8)  Population  size  and  system  capacity  are  infinite. 

(9)  The  system  does  not  allow  bulk  arrivals,  bulk  services,  balking  or 
reneging . 

(10)  The  system  operates  under  the  steady-state  conditions  which  are 
assumed  to  exist  and  which  will  be  determined  in  the  sequel. 

1 . 4 Objectives  of  Research 

Our  model  is  similar  to  the  one  outlined  by  Bell[5],  but  the 

following  refinements  point  out  the  differences  between  our  work  and  the 

aforementioned  work  including  Bell's,  thus  justifying  our  research. 

(1)  Rather  than  be  concerned  with  deriving  the  optimal  operating 

policy  structure  explicitly,  one  of  the  objectives  of  this  research 
is  to  derive  the  system  characteristics  which  might  be  incorporated 
in  a total  expected  cost  function  without  using  the  explicit  form 
of  the  steady- state  joint  probability  mass  function  of  the  number 
of  customers  and  the  number  of  operating  service  stations  in  the 
system.  Our  approach  could  be  applicable  to  find  out  the  system 
characteristics  of  other  Markovian  queueing  problems  in  which  there 
are  some  difficulties  in  obtaining  explicitly  the  steady-state 
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joint  probability  mass  function  of  the  number  of  customers  and  the 
number  of  operating  service  stations  in  the  system. 

(2)  Another  objective  of  this  research  is  to  derive  important  system 
characteristics  which  might  not  be  derivable  from  the  steady- state 
joint  probability  mass  function  of  the  number  of  customers  and  the 
number  of  operating  service  stations  in  the  system.  This  may 
involve  for  example  the  determination  of  the  probability  mass 
function  of  the  number  of  activations  or  removals  of  service 
stations  during  a busy  period,  the  Laplace  transforms  of  the 
probability  density  functions  of  the  busy  period,  the  length  of 
time  when  one  service  station  is  operating  and  the  length  of  time 
when  both  service  stations  are  operating  simultaneously  during  the 
busy  period.  These  system  characteristics  might  be  essential  to 
provide  fundamental  information  for  analysis  of  the  system 
behavior.  In  particular,  a number  of  existing  results,  namely,  the 
Laplace  transforms  of  the  probability  density  functions  of  the  busy 
periods  for  the  ordinary  M/M/2  system  derived  by  Conolly[8],  the 
M/M/1  system  operating  under  the  N policy  originated  by  Yadin  and 
Naor[27]  and  the  ordinary  M/M/1  system  derived  by  Conolly[8]  or 
Kleinrock[15]  can  be  recovered  from  the  obtained  Laplace  transform 
of  the  probability  density  function  of  the  busy  period  for  our 
model . 

(3)  Even  though  the  main  purpose  of  removable  service  stations  in 
queueing  system  is  to  utilize  the  idle  service  stations,  most  of 
the  total  expected  cost  functions  developed  in  the  area  of  multi- 
removable service  stations  are  not  concerned  with  the  utilization 
of  the  idle  service  stations  as  a possible  source  of  cost  element. 
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Furthermore,  the  same  cost  structures  for  each  operating  service 
station  are  assumed  regardless  of  changes  in  the  number  of 
operating  service  stations.  Because  of  these  assumptions,  lack  of 
applicability  for  the  total  expected  cost  function  to  various 
situations  is  unavoidable.  Hence,  total  expected  cost  function  per 
unit  time,  which  is  another  objective  of  this  research,  is 
constructed  by  considering  the  utilization  of  the  idle  service 
stations  as  a possible  source  of  cost  element  and  assuming 
different  structures  for  several  cost  elements  based  on  the  number 
of  operating  service  stations.  Thus,  this  total  expected  cost 
function  per  unit  time  could  be  applicable  to  a variety  of 
situations . 


1 . 5 Overview  of  the  Next  Chapters 

In  Chapter  2,  we  derive  the  system  characteristics  which  might  be 
incorporated  in  a total  expected  cost  function  using  the  explicit  form 
of  the  steady- state  joint  probability  density  function  of  the  number  of 
customers  and  the  number  of  operating  service  stations  in  the  system. 

The  advantage  of  this  approach  is  to  obtain  certain  system 
characteristics  relatively  easily  from  well  known  formulas.  However, 
although  the  explicit  form  of  this  joint  probability  mass  function  is 
known,  nevertheless  it  is  very  difficult  to  find  from  it  the  probability 
mass  function  of  the  number  of  activations  or  removals  of  service 
stations  during  the  busy  period,  the  probability  density  function  of  the 
busy  period  and  so  on. 

In  order  to  overcome  the  weakness  of  the  approach  used  in  Chapter 
2,  we  analyze  in  Chapter  3 the  busy  period  to  provide  more  important 
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information  of  the  system  behavior.  The  results  of  the  analysis  of  the 
busy  period  provide  the  fundamental  relations  to  obtain  the  probability 
mass  function  of  the  number  of  activations  or  removals  of  service 
stations  during  the  busy  period,  the  Laplace  transform  of  the 
probability  density  function  of  the  busy  period  and  so  on. 

In  Chapter  4,  we  derive  the  probability  mass  function  of  the  number 
of  activations  or  removals  of  service  stations  during  the  busy  period 
using  the  results  of  the  analysis  of  the  busy  period  in  Chapter  3 and 
the  properties  of  the  classical  gambler's  ruin  problem.  The  explicit 
form  of  this  probability  mass  function  plays  key  role  to  proceed  towards 
other  objectives  of  this  research. 

In  Chapter  5,  we  derive  the  Laplace  transforms  of  the  probability 
density  functions  of  the  busy  period,  the  busy  cycle,  the  length  of  time 
when  one  service  station  is  operating  and  the  length  of  time  when  both 
service  stations  are  operating  simultaneously  during  a busy  period.  We 
also  use  the  properties  of  the  classical  gambler's  ruin  problem  along 
with  the  results  of  Chapter  2 and  3.  We  recover,  as  a special  case,  the 
Laplace  transforms  of  the  probability  density  functions  of  the  busy 
period  for  the  ordinary  M/M/2  system,  the  M/M/1  system  operating  under 
the  N policy  and  the  ordinary  M/M/1  system. 

In  Chapter  6,  we  develop  alternative  methodologies  to  find  certain 
system  characteristics  that  will  be  used  to  construct  a total  expected 
cost  function  in  Chapter  7 based  on  the  results  of  Chapter  3,  4 and  5. 

In  Chapter  7,  we  construct  a total  expected  cost  function  per  unit 
time  based  on  the  results  of  the  previous  chapters . We  assume  that  the 
cost  elements  include  the  cost  for  the  service  stations  removed  from  the 
system,  the  holding  cost  for  waiting  customers,  the  operating  cost  for 
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operating  service  stations  and  the  switching  cost  for  activations  and 
removals  of  service  stations.  We  also  assume  (i)  different  operating 
cost  for  each  operating  service  stations  depending  on  the  number  of 
operating  service  stations  similar  to  the  concepts  of  individual  and 
joint  ordering  policies  in  inventory  theory  and  (ii)  the  switching  cost 
for  activations  and  removals  of  service  stations  is  not  necessarily 
incurred  whenever  service  stations  are  activated  or  removed.  The 
optimal  customer  levels  (K  at  which  service  stations  should  be 

activated  and  removed  are  determined  numerically  in  order  to  minimize 
the  total  expected  cost  per  unit  time.  We  also  perform  a sensitivity 
analysis  for  changes  of  the  optimal  customer  levels  along  with  changes 
in  specific  cost  elements. 

Finally,  in  Chapter  8,  we  discuss  extension  of  our  methodologies  to 
analyze  the  busy  period  for  other  Markovian  queueing  models.  Topics  for 
future  research  in  the  area  of  controllable  M/M/2  queueing  system  are 
also  suggested. 


CHAPTER  2 


DERIVATIONS  OF  SYSTEM  CHARACTERISTICS  USING  THE  JOINT 
PROBABILITY  MASS  FUNCTION  OF  THE  NUMBER  OF  CUSTOMERS  AND 
THE  NUMBER  OF  OPERATING  SERVICE  STATIONS  IN  THE  SYSTEM 

2 . 1 Introduction 

The  main  purpose  of  this  chapter  is  to  derive  the  important  system 
characteristics  for  a steady-state  M/M/2  queueing  system  with  removable 
service  stations  operating  under  the  (0,K,N,M)  policy  over  an  infinite 
time  horizon.  The  methodology  used  in  this  chapter  requires  the 
explicit  form  of  the  steady- state  joint  probability  mass  function  of  the 
number  of  customers  and  number  of  operating  service  stations  in  the 
system. 

The  following  six  system  characteristics  which  might  be 
incorporated  in  a total  expected  cost  function  per  unit  time  are  derived 
explicitly  in  this  chapter. 

(1)  Steady-State  Probability  that  ' i'  Service  Stations  are  Operating 

and  '1'  Customers  are  in  the  System,  denoted  by  P(i,i)  where  i»0,1.2  and 
■1  ~0 . 1 . 2 • • • : The  P(i,j)'s  are  obtained  by  solving  a set  of  steady-state 

difference  equations  recursively. 

(2)  Probability  Generating  Function  of  the  Number  of  Customers  in  the 

System,  denoted  by  G(z)  : G(z)  is  derived  explicitly  using  the  set  of 

steady-state  difference  equations. 
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(3)  Expected  Number  of  Customers  in  the  System,  denoted  by  EfJ]  : 
From  the  explicit  form  of  G(z) , E[J]  is  derived  by  using  the  property 

dG(z) 


E[J] 


dz 


z-1 . 


(^)  Expected  Number  of  Operating  Service  Stations,  denoted  by  EfOl  : 
Using  the  well-known  Little's  formula,  E[0]  is  obtained. 

(5)  Expected  Length  of  the  Idle  Period,  the  Busy  Period  and  the  Busy 
Cycle,  denoted  by  E[I],  E[B1  and  EfCT,  respectively  : Define  the  idle 

period,  the  busy  period  and  the  busy  cycle  as  following. 

1)  Idle  period  denoted  by  I ; the  length  of  time  when  both  service 
stations  remain  inoperative  temporarily. 

2)  Busy  period  denoted  by  B ; the  length  of  time  when  at  least  one 
service  station  is  operating. 

3)  Busy  cycle  denoted  by  C ; the  length  of  time  from  the  beginning 
of  the  last  idle  period  to  the  beginning  of  the  next  idle 
period.  In  other  words,  the  busy  cycle  is  the  sum  of  the  idle 
period  and  the  busy  period. 

When  both  service  stations  remain  inoperative  temporarily,  no  service 
stations  are  allowed  to  operate  until  the  number  of  customers  who  are 


waiting  for  service  in  the  system  reaches  N.  Hence,  the  idle  period  can 
be  represented  by  the  sum  of  the  Aj's  where  j=l,2,---,N  which  are  the 
(j-l)th  and  the  j th  customers'  interarrival  times.  By  the  assumptions 
of  system  operation,  the  A j ' s are  identically  and  independently 
distributed  with  E[Aj]  - 1/A.  Therefore, 

E[I]  - N E[A.]. 
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Let  P[0-0]  be  the  probability  that  no  service  stations  are  operating. 
Then,  P[0— 0]  can  be  computed  from 

N-1 

P[0-0]  - S P(0,j). 
j“0 

Since 

E[I] 

P[0-0] , 

E[I]  + E[B] 

then,  E[B]  is  obtained  from 

E[I]  - E[I]  P[0-0] 

E[B] . 

P[0-0] 

Thus,  E[C]  can  be  derived  explicitly  from 

E[C]  - E[I]  + E[B]  . 

(6)  Expected  Number  of  Activations  (Removals)  of  Service  Stations  per 
Time,  denoted  by  E[R^]  (E[R.^])  : Define  the  random  variables 

and  as  the  number  of  activations  and  removals  of  service  stations  per 
unit  time,  respectively.  In  any  busy  cycle,  one  of  the  service  stations 
removed  from  the  system  should  be  activated  immediately  whenever  a new 
customer  arrives  at  the  system  (i)  while  N-1  customers  are  in  the  system 
and  both  service  stations  remain  inoperative  temporarily,  or  (ii)  while 
M-1  customers  are  in  the  system  and  one  service  station  is  operating. 
Since  the  assumed  customers'  mean  arrival  rate  is  A,  E[R  1 is  obtained 
from 

E[R^]  - AP(0,N-1)  + AP(1,M-1). 

Similarly,  E[Rj_]  is  obtained  form 


E[R^]  = mP(1,1)  + 2mP(2,K+1). 
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2 . 2 Steady- State  Probability  that  'i'  Service  Stations  are  Operating 
and  ' i ' Customers  are  in  the  System 


Since  we  have  assumed  that  the  customers'  interarrival  times  and 
the  service  times  of  two  undistinguishable  service  stations  are 
independently  and  exponentially  distributed,  the  number  of  customers  in 
the  system  has  Markovian  properties.  From  these  Markovian  properties, 
we  can  deduce  that  the  mean  transition  rate  into  a given  state  is  equal 
to  the  mean  transition  rate  out  of  that  state.  Hence,  we  can  write  a 
set  of  steady- state  probability  equations  which  will  lead  to  expressions 
for  the  probabilities  of  being  in  any  of  admissible  states.  Define 
P(i.j)  where  i-0,1,2.  and  J-0 , 1 , 2 , 3 , • • • as  the  steady- state  probability 
that  'i'  service  stations  are  operating  and  'j'  customers  are  in  the 
system.  Referring  to  the  rate  diagram  for  the  M/M/2  system  operating 
under  the  (0,K,N,M)  policy  appearing  in  Figure  2.1,  the  admissible 
states  of  the  system  can  be  partitioned  into  several  disjoint  subsets 


which  lead  to  the  following  equations: 

AP(i,j)  - AP(i.j+l)  (i-0,  0<j<N-2),  (2.1) 

AP(i,j)  = /iP(i+l,j-i-l)  (i-0,  j-0),  (2.2) 

(A+/i)P(i,j)  = MP(i,j+l)  (i-1,  j-1),  (2.3) 

(A+/i)P(i,j)  - /iP(i,j-i-l)  + AP(i,j-l)  (i=l,  2<j^-2,  j/K,N),  (2.4) 

(A+;i)P(i,j)  - ^P(i,j+1)  + AP(i,j-l)  + 2^P(i-t-l, j+1)  (i-1,  j-K),  (2.5) 

(A-t-/i)P(i,j)  - n?(i,2+l)  + AP(i,j-l)  + AP(i-l,j-l)  (i-1,  j-N),  (2.6) 

(A+^)P(i,j)  = AP(i,j-l)  (i-1,  j-M-1),  (2.7) 

(A+2p)P(i,j)  - 2AtP(i,j+l)  (i-2,  j=K+l)  , (2.8) 

(A-t-2^)P(i,  j ) - 2AtP(i,  j+1)  + AP(i,j-l)  (i=2,  K+2<j  , j/M),  (2.9) 


(A+2/i)P(i,j)  = 2/iP(i,j+l)  + AP(i,j-l)  -i-  AP(i-l,j-l)  (i=2,j=M).  (2.10) 
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Figure 


number  of  operating  service  stations 


2 10 

J I L 


2.1  The  rate  diagram  of  the  model  with  an  infinite 
system  capacity  and  an  infinite  population  size 
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Define  p - \/p.  P(i,j)  for  i-0,1,2.  and  j=l , 2 , 3 , • • • are  obtained 
in  terms  of  P(0,0)  by  solving  equations  (2.1)  to  (2.10)  recursively. 
The  summary  of  the  results  is  given  below. 

For  p 1 

P(0,0)  (i-0,  l<j<N-l) 

pa-ph 

P(0,0)  (i-1,  l<j^) 

(1-p) 


P(iJ) 


p(l-p^)  p(p^-p“-N)(l-pJ-K) 

+ ' 

(1-p)  (l-p)d-p^-K) 


P(0,0)  (i-1,  K+l<j:^) 


p(l-p^)(pJ-N-p«-N) 

— — P(0,0)  (i-1,  N+l<j<M-l)  (2.11) 

(l-p)d-p^-K) 

pM-N+l(i.^N)  j;L-(p/2)j-^] 

— P(0,0)  (i-2,  K+l<j^) 

(2-p)(l-pM-^<) 

P^'N+1(1-/)  [1-  (p/2)^-‘^]  (p/2)j 

— P(0,0)  (i-2,  j^+1), 

(2-p)d-pM-i^) 


and  for  p - 1 


P(i.J) 


P(0,0)  (i-0,  l<j<N-l) 

jP(0,0)  (i-1,  l<j<K) 

NK  + (M-K-N)j 

P(0,0)  (i-1,  K+1<J^) 

M-K 

N(M-j) 

P(0,0)  (i-1,  N+l<j<M-l) 

M-K 

N[l-(l/2)j-^] 

P(0,0)  (i-2,  K+l<j^) 

M-K 


(2.12) 
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P(i.j)  - 


M-K 


P(0,0) 


(i=2,  j^+1). 


(2.12) 


Evaluation  of  P(0,0)  can  be  done  by  the  normalizing  condition 

N-1  M-1  00 

S P(0,j)  + E P(l,j)  + S P(2,j)  - 1. 

j-0  j-1  j-K+1 

In  order  for  the  left  hand  side  of  equation  above  to  converge,  p should 
be  less  than  two.  The  normalizing  condition  can  also  be  obtained  from 
the  probability  generating  function  of  the  number  of  customers  in  the 
system.  The  evaluation  of  P(0,0)  is  deferred  until  the  probability 
generating  function  of  the  number  of  customers  in  the  system  is  derived 
in  terms  of  P(0,0). 


2 . 3 Probability  Generating  Function  of 

the  Number  of  Customers  in  the  System 

Define  the  probability  generating  function  of  the  number  of 

customers  in  the  system,  G(z)  as 

N-1  . M-1  . 00 

G(z)  - E zJ  P(0,j)  + E zj  P(l,j)  + E zJ  P(2,j).  (2.13) 

j-0  j-1  j=K+l 

Using  the  steady-state  probability  equation  (2.1),  multiply  by  z^ 

and  sum  over  all  possible  ^values  of  j . Then, 

N-1  . N-1 

E zjp(O.j)  - E zjp(0,0) 
j-0  j-0 


1-z 


P(0,0). 


(2.14) 
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Similarly,  using  equations  (2.2)  to  (2.7), 


M 


pzP(0,0)  + E zJ(l+^)P(l,j-l) 

j-2 


M-1  M 

+ S zJp(l,J)  + S pzjp(l,N-2) 

j-1  j-3 

+ 2z*^'^^P(2,K+l)  + pz^'^^P(0,N-l) . 


Hence , 
M-1 


S zJp(l,j)  - 

J-1 


pzP(0,0)  - pz^‘^^P(0,0) 
(l-pz)d-z) 


pz^‘^^P(l,M-l)  - 2z^'^^P(2,K+l) 
(l-pz)d-z) 


(2.15) 


Finally,  for  S z3p(2,j),  using  equations  (2.8)  to  (2.10) 
j-K+1 

yields 

CO  00 

S zj(l+p/2)P(2,j)  - S zip(2,j+l) 

j=K+l  j=K+l 


+ S zi(p/2)P(2,j-l)  + z^(p/2)P(l,M-l) 
j-K+2 


Thus , 


S zJ  P(2,j)  - 
j-K+1 


2z*^'^^P(2,K+l)  - pz^'^^  P(1,M-1) 


(2-pz)(l-z) 


(2.16) 


Substituting  (2.14),  (2.15)  and  (2.16)  into  relation  (2.13)  and 
simplifying  yield 


(1-zN) 


G(z) 


(l-z)d-pz) 


P(0,0)  + 


pz^+^P(l,M-l)  - 2z^'^^P(2,K+l) 
(l-z)(l-pz)(2-pz) 


(2.17) 
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From  the  steady- state  probabilities  obtained  in  (2.11),  P(1,M-1)  and 
P(2,K+1)  are  given  by 


P(l.M-l) 


P(2,K+1)  - 


P(0,0), 


P(O.O) 


(2.18) 


(2.19) 


Substituting  the  results  for  P(1,M-1)  and  P(2,K+1)  in  (2.18)  and  (2,19) 

into  (2.17)  and  simplifying  yield 

[(l-zN)(2-pz)(l-p^-K)  + p^-N+1(1-/)(zM+1-2K+1)]p(o,0) 

G(z) — — . (2.20) 

(l-p'^-*^)(l-z)(l-pz)(2-pz) 

The  nonnalizing  condition  is  obtained  by  using  the  property 

- 1. 


G(z) 

Apply  L'Hospital's  rule  twice, 
G(z) 


z-1 


-N(l-p^-^)(2-p)  + p^-N+l(l-pN)(f^.K) 


z=l 


Hence , 


P(0,0) 


-d-p^*^)(l-p)(2-p) 


(l-p^-^)(l-p)(2-p) 


P(0,0) 


N(l-p^-^)(2-p)  - p^-N+l(l-pN)(M.K) 


for  p / 1 


(2.21) 


N(M+K-N+4) 


for  p - 1. 


By  substituting  the  value  above  for  P(0,0)  in  (2.21)  into  equation 
(2.20),  it  immediately  follows  that 


G(z) 


(l-p)(2-p)[(l-pM-K)(l-zN)(2-pz)  - p^'N'^l(l-pN)(zM+1.2K+l)j 
[N(l-p^'^)(2-p)  - p^'^"^^(l-p^)(M-K)] (l-z)(l-pz)(2-pz) 


(2.22) 
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2 . 4 Expected  Number  of  Customers  in  the  System 


Let  J be  the  number  of  customers  in  the  system.  E[J]  can  be 
obtained  explicitly  from 


dG(z) 

E[J] 

dz 


z-1 . 


Hence,  we  obtain  E[J]  as 


for  p / 1 

M+K+1  p(3-2p)  N(1-p^‘^){(M+K-N)(2-p)+4) 

E[J] + ^ 

2 (l-p)(2-p)  2{N(l-p^*^)(2-p)  - p^*N+l(M-K)(l-pN) 

(2.23) 

for  p - 1 

M+3  6(N-l)+(K+N)(K-N)+M(N-4) 

E[J] + (2.24) 

3 3(M+K-N+4) 


2 . 5 Expected  Number  of  Operating  Service  Stations 


Under  steady-state  conditions,  define  the  random  variables  0 J W 

’ n * 


and  Wq  as 


1)  0 ; the  number  of  operating  service  stations, 

2)  : the  number  of  customers  who  are  waiting  for  service  in  the 

system  ( does  not  include  the  customers  who  are  being 
served. ) , 

3)  W : the  customer's  waiting  time  in  the  system  ( W does  include 

service  time.), 

4)  Wq  ; the  customer's  waiting  time  for  service  ( does  not 

include  the  service  time.). 
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Then,  the  following  relations  hold: 

E[J]  - E[Jq]  + E[0] . 

E[W]  - E[Wq]  + E[SJ. 

From  Little's  formulas,  we  have 

E[J]  - AE[W]  and  E[Jq]  - AE[W^] . 

Hence,  we  obtain  E[0]  as 

E[0]  - E[J]  - E[Jq] 

- A{E[W]  - E[Wq]} 

- AE[Sj^] 

- P.  (2.25) 

Even  though  E[J^] , E[W]  and  E[W^]  are  not  incorporated  in  the  total 
expected  cost  function  per  unit  time,  they  can.be  obtained  explicitly 
from  the  above  relations  and  the  Little's  formulas  such  as 


E[Jq]  - 

E[J]  - E[0], 
1 

(2.26) 

E[W]  - 

E[J] , 

A 

1 

(2.27) 

E[Wq]  - 

— E[J^]. 

(2.28) 

2 . 6 Expected  Length  of  the  Idle  Period, 
the  Busy  Period  and  the  Busy  Cycle 

This  removable  queueing  system  alternates  between  durations  when 
both  service  stations  remain  inoperative  temporarily  and  durations  when 
at  least  one  service  station  is  operating.  To  clarify  the  system 
states,  the  idle  period,  the  busy  period  and  the  busy  cycle  are  defined 


as  follows. 
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1)  Idle  period  denoted  by  I : the  length  of  time  when  both  service 

stations  remain  inoperative  temporarily;  this  begins  the  instant  when 
the  last  customer  departs  from  the  system  and  endures  until  the  arrival 
of  the  next  Nth  customer. 

2)  Busy  period  denoted  by  B : the  length  of  time  when  at  least  one 

service  service  station  is  operating;  this  begins  the  instant  when  the 
first  service  station  starts  operating  at  which  time  the  number  of 
waiting  customers  reaches  N and  both  service  stations  remain  inoperative 
temporarily,  and  lasts  until  there  are  no  customers  in  the  system.  In 
other  words,  the  busy  period  is  the  length  of  time  between  the  end  of 
the  last  idle  period  and  the  beginning  of  the  next  idle  period. 

3)  Busy  cycle  denoted  by  C ; the  length  of  time  between  the  beginning 
of  last  idle  period  and  the  beginning  of  the  next  idle  period.  In  other 
word,  the  busy  cycle  is  the  sum  of  the  idle  period  and  the  busy  period 
such  as 

C = I + B,  (2.29) 

or 

E[C]  - E[I]  + E[B] . (2.30) 

After  both  service  stations  become  inoperative  temporarily,  no 
service  stations  may  operate  until  the  number  of  customers  who  are 
waiting  for  service  in  the  system  reaches  N.  Hence,  the  idle  period  can 
be  represented  by  the  sum  of  the  A j ' s (j-l,2---,N)  which  are  the  (j-l)th 
and  the  jth  customers'  interarrival  times;  thus, 

I = + A2  + A3  + • • • + Ajj_3  + Ajj. 


(2.31) 
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number  of 
customers 


Figure  2.2  The  sample  function  of  the  busy  cycle 


By  the  basic  assumptions  of  system  operation,  the  A j ' s are  identically, 
independently  and  exponentially  distributed  random  variables  with 
E[Aj]  - 1/A.  Thus,  the  expected  length  of  idle  period,  E[I]  is  obtained 
from 


N 

E[I]  - S E[A^]  = N E[A.]  - N/A.  (2.32) 

j=l  J J 

Let  P[0=0]  be  the  probability  that  no  service  stations  are 
operating.  Then,  P[0—0]  can  be  computed  from 

N-1 

P[0-0]  - S P(0,j). 
j-0 

Hence , 


P[0-0] 


N(l-p)(2-p)(l-p^-^) 
N(l-p”-^)(2-p)  - p^'N+l(l-pN)(f^.K) 
2 


for  p / 1 

(2.33) 


M+K-N+4 


for  p - 1. 


36 


Since 

E[I]  E[I] 

P[0=0]  , 

E[C]  E[I]  + E[B] 

then 

E[I]  - E[I]  P[0-0] 

E[B] . 

P[0=0] 

Therefore , 


for  p \ 

(2.34) 


1 

A 


E[B] 


Np 


(1-p)  (l-p)(2-p)(l-p”-‘^) 


N(M+K-N+2) 


2A 


for  p = 1. 


Substituting  the  explicit  forms  for  E[I]  in  (2.32)  and  E[B]  in 
(2.34)  into  relation  (2.30)  yields 


E[C] 


N p^'^'^^(l-p^)  (M-K) 

(1-p)  (l-p)(2-p)(l-p^-^) 

N(M+K-N+4) 

for  p = 1. 

2A 


1 

A 


for  p / 1 


(2.35) 


2 . 7 Expected  Number  of  Activations  and  Removals  of 
Service  Stations  per  Unit  Time 


Define  the  random  variables  and  R^.  as  positive  integers  to  be 
the  number  of  activations  and  removals  of  service  stations  per  unit 


time . 
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One  of  the  service  stations  removed  from  the  system  should  be 
activated  immediately  (i)  whenever  the  number  of  customers  in  the  system 
reaches  N while  both  service  stations  remain  inoperative  temporarily,  or 
(ii)  whenever  the  number  of  customers  in  the  system  reaches  M while  one 
service  stations  are  operating  and  the  other  remains  inoperative.  Since 
the  assumed  customers'  mean  arrival  rate  per  is  A,  E[R_]  is  obtained 
from 


E[R^]  = AP(0,N-1)  + AP(1,M-1).  (2.37) 

Substituting  P(0,N-1)  and  P(1,M-1)  in  (2.11)  for  / 1 and  in  (2.12) 
for  p - 1 with  P(0,0)  in  (2.21)  and  simplifying  yield 


A(l-p)(2-p)[(p^-N-pM)  + (1-pM-K)] 

N(2-p)(l-p‘*^-^)  - 


for  p / 1 


(2.38) 


2A(M+N-K) 


N(M-K) (M+K-N+4) 


for  p = 1. 


Similarly,  the  service  station  just  finishing  service  should  be 
removed  from  the  system  immediately  (i)  whenever  the  number  of  customers 
in  the  system  reaches  0 while  one  service  station  is  operating  and  the 
other  is  removed  from  the  system,  or  (ii)  whenever  the  number  of 
customers  in  the  system  reaches  K while  both  service  stations  are 
operating  simultaneously.  However,  mean  service  rate  of  each  operating 
service  station  per  unit  time  is  p.  Thus,  E[Rj_]  can  be  obtained  from 

E[R^]  = mP(1,1)  + 2pP(2,K+l).  (2.39) 

Substituting  the  explicit  values  of  P(l,l)  and  P(2,K+1)  in  (2.11)  and 
(2.12)  into  relation  (2.39)  yields  E[Rj_]  explicitly  which  is  the  same  as 
the  results  of  E[R  ] obtained  in  (2.38). 
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In  the  next  chapter,  we  will  analyze  the  busy  period  systematically 
to  provide  basic  information  to  derive  the  probability  mass  function  of 
the  number  of  activations  or  removals  of  service  stations  during  the 
busy  period,  the  Laplace  transform  of  the  probability  density  function 
of  the  busy  period  and  so  on. 


CHAPTER  3 


ANALYSIS  OF  THE  BUSY  PERIOD 
3 . 1 Introduction 

In  this  chapter,  we  analyze  the  busy  period  systematically  to 
provide  basic  information  for  (1)  obtaining  the  probability  mass 
function  of  the  number  of  activations  or  removals  of  service  station 
during  the  busy  period,  (2)  deriving  the  Laplace  transform  of  the 
probability  density  function  of  the  busy  period,  (3)  deriving  the 
Laplace  transforms  of  the  probability  density  functions  of  the  length  of 
time  when  one  service  station  is  operating  and  of  the  length  of  time 
when  both  service  stations  are  operating  simultaneously  and  (4) 
developing  alternative  methodologies  to  obtain  certain  system 
characteristics  that  will  be  used  to  construct  a total  expected  cost 
function  in  Chapter  7,  without  using  the  explicit  form  of  P(i,j). 

The  methodology  used  in  this  chapter  is  to  partition  the  busy 
period  into  several  random  intervals  depending  on  the  numbers  of 
customers  and  the  number  of  operating  service  stations  in  the  system. 

In  other  words,  any  busy  period  can  be  represented  in  terms  of  these 
random  intervals.  Derivation  of  the  relation  between  the  busy  period 
and  the  random  intervals  is  a prime  objective  of  this  chapter.  We  also 
analyze  each  random  interval  based  on  the  fluctuation  of  the  number  of 
customers  in  the  system.  The  properties  of  the  random  walk  or  the 
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classical  gambler's  ruin  problem  are  used  to  characterize  the 
fluctuation  of  the  number  of  customers  in  the  system. 

3 . 2 Preliminary  Work  for  Analysis  of  the  Busy  Period 

To  proceed  towards  a formal  analysis  of  the  busy  period,  define  the 
random  variables  as 

1)  and  F^.  ; both  positive  integers  to  be  the  number  of  activations 
and  removals  of  service  stations  during  the  busy  period, 

2)  D ; a nonnegative  integer  to  be  the  number  of  times  the  level  of 
customers  in  the  system  reaches  M while  one  service  station  is 
operating  during  the  busy  period. 

The  main  objective  of  this  section  is  to  show  that 

i)  F^  - F^. 

ii)  F^  = D + 1. 

In  any  busy  cycle,  termination  of  the  idle  period  and  initiation  of 
the  busy  period  occur  at  the  same  time  when  the  number  of  customers  who 
are  waiting  for  service  reaches  N while  both  service  stations  remain 
inoperative  temporarily.  This  indicates  that  the  busy  period  is 
initiated  when  one  of  the  two  service  stations  removed  from  the  system 
starts  operating  and  lasts  until  both  service  stations  become 
inoperative  temporarily  again.  Since  the  system  operates  under  steady- 
state  conditions,  any  service  station  which  is  operating  during  the  busy 
period  should  be  removed  before  the  beginning  of  the  next  idle  period. 
Thus,  the  number  of  activations  of  service  stations  is  equal  to  the 
number  of  removals  of  service  stations  during  the  busy  period.  Hence, 


we  have 
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Fa  " Fr-  (3.1) 

When  one  service  station  is  operating  and  one  service  station  is 
removed  from  the  system,  the  level  of  customers  may  fluctuate  between  1 
and  (M-1)  without  having  the  second  service  station  activated,  or  after 
reaching  M,  the  level  fluctuates  until  it  reaches  K at  which  time  only 
one  service  station  is  operating,  after  which  time  the  level  M is 
reached  again  and  so  on.  Hence,  there  are  no  specific  upper 
restrictions  of  fluctuations  in  the  level  of  customers  when  both  service 
stations  are  operating  simultaneously.  This  indicates  that  the  number 
of  customers  in  the  system  can  be  allowed  to  increase  indefinitely 
without  changing  the  number  of  operating  service  stations.  Therefore, 
the  activations  or  the  removals  of  the  service  stations  depend  on  when 
the  number  of  customers  in  the  system  reaches  specified  levels  and  how 
many  service  stations  are  operating.  Thus,  the  following  theorem  might 
be  stated  based  on  the  fluctuation  of  the  number  of  customers  in  the 
system  depending  on  the  number  of  operating  service  stations. 

Theorem  3.1 

The  number  of  activations  or  removals  of  the  service  stations 
during  the  busy  period  is  equal  to  D + 1,  i.e., 

- D + 1 for  all  D.  (3.2) 

Proof  : Depending  on  the  values  of  D,  one  of  the  following  cases  can 

occur  during  the  busy  period. 

Case  1 : D — 0 

Only  the  service  station  which  is  activated  at  the  beginning  of  the  busy 
period  is  operating  during  the  entire  busy  period.  The  other  available 
service  station  remains  removed  from  the  system  during  the  entire  busy 
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period.  This  case  occurs  only  when  the  number  of  customers  in  the 
system  is  always  less  than  M during  the  busy  period.  In  this  case,  the 
number  of  activations  or  removals  of  the  service  stations  during  the 
busy  period  is  1,  i.e., 


Case  2 : D - 1 

One  of  the  two  service  stations  removed  from  the  system  starts  operating 
at  the  beginning  of  the  busy  period.  The  remaining  service  station 
removed  from  the  system  also  starts  operating  when  the  number  of 
customers  in  the  system  reaches  M during  the  busy  period.  Hence,  the 
number  of  operating  service  stations  becomes  two.  At  a later  time,  when 
the  number  of  customers  in  the  system  reaches  K,  the  service  station 
just  finished  service  is  removed  from  the  system  immediately.  After 
this  time,  the  number  of  customers  in  the  system  never  reaches  M again. 
Therefore,  the  busy  period  is  terminated  when  the  remaining  service 
station  is  removed  from  the  system.  In  this  case,  the  number  of 
activations  or  removals  of  service  stations  during  the  busy  period  is  2, 
i.e.  , 


Case  3 : D > 2 

One  of  the  two  service  stations  removed  from  the  system  starts  operating 
at  the  beginning  of  the  busy  period.  The  remaining  service  station 
removed  from  the  system  also  starts  operating  when  the  number  of 
customers  in  the  system  reaches  M during  the  busy  period.  Hence,  the 
number  of  operating  service  stations  becomes  two  for  the  first  time 
during  the  busy  period.  At  a later  time,  when  the  number  of  customers 
in  the  system  reaches  K,  the  service  station  just  finished  service  is 
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removed  immediately.  Thus,  the  number  of  operating  service  stations 
becomes  one  again.  In  other  words,  when  the  number  of  customers  in  the 
system  reaches  M from  K without  reaching  0 while  one  service  station  is 
operating,  the  number  of  operating  service  stations  becomes  two.  Then, 
when  the  number  of  customers  in  the  system  reaches  K from  M without 
reaching  0,  the  number  of  operating  service  stations  becomes  one  again. 
Thus,  the  number  of  operating  service  stations  alternates  between  one 
and  two  (D-1)  times.  After  that,  the  number  of  customers  in  the  system 
never  reaches  M again.  Hence,  the  busy  period  is  terminated  when  the 
remaining  operating  service  station  is  removed.  In  this  case,  the 
number  of  activations  or  removals  of  the  service  stations  during  the 
busy  period  is  (D  + 1) , i.e., 

Fr  - Fa  - D + 1. 

Hence,  from  the  above  three  possible  cases,  the  theorem  follows 
immediately.  Q.E.D. 

From  (3.2),  we  have 

E[Fj.]  - E[F^]  - E[D]  + 1.  (3.3) 

Next,  we  will  partition  the  busy  period  into  the  random  intervals 
based  on  the  number  of  customers  and  the  number  of  operating  service 
stations  in  the  system. 

3 . 3 Partition  of  the  Busy  Period 


In  order  to  represent  mathematically  the  fluctuations  of  the  number 
of  customers  in  the  system  during  the  busy  period,  define  the  following 


random  time  intervals. 
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1)  Uq  : the  length  of  time  from  the  beginning  of  the  busy  period  to 

the  first  time  the  number  of  customers  in  the  system  reaches  0 without 
reaching  M.  Hence,  Uq  is  the  first  passage  time  from  N to  0 without 
reaching  M.  During  Uq,  the  level  of  customers  in  the  system  may 
fluctuate  between  1 and  (M-1)  before  it  reaches  0 without  reaching  M 
given  that  one  service  station  is  operating. 

2)  Uj^  : the  length  of  time  from  the  beginning  of  the  busy  period  to 

the  first  time  the  number  of  customers  in  the  system  reaches  M without 
reaching  0.  Hence,  U^  is  the  first  passage  time  from  N to  M without 
reaching  0.  During  U^,  the  level  of  customers  in  the  system  may 
fluctuate  between  1 and  (M-1)  before  it  reaches  M without  reaching  0 
given  that  one  service  station  is  operating. 

3)  U2  : the  length  of  time  from  the  instant  when  the  number  of 

customers  in  the  system  reaches  K to  the  end  of  the  busy  period  without 
reaching  M.  Hence,  U2  is  the  first  passage  time  from  K to  0 without 
reaching  M.  During  U2,  the  level  of  customers  in  the  system  may 
fluctuate  between  1 and  (M-1)  before  it  reaches  0 without  reaching  M 
given  that  one  service  station  is  operating. 

where  d-1 , 2 , • • • ,D-1  and  D > 1 : the  dth  length  of  time  from  the 
instant  when  the  number  of  customers  in  the  system  reaches  K while  both 
service  stations  are  operating  to  the  first  time  the  number  of  customers 
in  the  system  reaches  M without  reaching  0.  Hence,  is  the  first 
passage  time  from  K to  M without  reaching  0.  During  the  Y^j's,  the  level 
of  customers  in  the  system  may  fluctuate  between  1 and  (M-1)  before  it 
reaches  M without  reaching  0 given  that  one  service  station  is 
operating.  Note  that  the  Y^'s  are  identically  and  independently 
distributed  and  independent  of  D and  F . 
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Figure  3.1  The  sample  functions  of  the  fluctuations  of 
the  number  of  customers  in  the  system 
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5)  where  d-l,2,---,D  and  D > 1 : the  dth  length  of  time  from  the 

Instant  when  the  number  of  customers  in  the  system  reaches  M while  one 
service  station  is  operating  to  the  first  time  the  number  of  customers 
in  the  system  reaches  K.  Hence,  is  the  first  passage  time  from  M to 
K.  During  the  V^'s,  the  level  of  customers  in  the  system  may  fluctuate 
between  (K+1)  and  an  indefinite  level  which  is  greater  than  or  equal  to 
M before  it  reaches  K given  that  both  service  stations  are  operating 
simultaneously.  Note  that  the  V^'s  are  identically  and  independently 
distributed  and  independent  of  D and  F^. 

We  note  in  particular  that  the  length  of  any  busy  period  may  be 
represented  in  terms  of  the  random  variables  Uq,  U^,  U2,  and  for 
all  d. 


Theorem  3 , 2 

Any  busy  period,  B can  be  expressed  in  the  following  form: 


B 


Uq  for  D - 0 

+ U2  for  D - 1 

D-1  D 

Uf  + S + S + U2  for  D > 2. 


(3.4) 


Proof  : Theorem  3.2  follows  the  interpretation  of  the  sample  functions 

of  the  fluctuations  of  the  number  of  customers  in  the  system  depending 
on  the  number  of  operating  service  stations  appearing  in  Figure  3.1. 


Theorem  3.2  is  the  basic  information  for  (i)  obtaining  the 
probability  mass  function  of  F^,  (ii)  deriving  the  Laplace  transform  of 
the  probability  density  function  of  the  busy  period,  B and  (iii) 
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developing  alternative  methodologies  to  obtain  certain  system 
characteristics  without  using  the  explicit  form  of  P(i,j). 

However,  during  Uq,  U^,  U2  and  for  all  d,  exactly  one  service 
station  is  operating.  During  for  all  d,  both  service  stations  are 
operating  simultaneously.  Hence,  any  busy  period  can  be  represented  by 
the  sum  of  the  length  of  time  when  one  service  station  is  operating  and 
the  length  of  time  when  both  service  stations  are  operating 
simultaneously.  Let  and  B2  be  the  random  variables  which  represent 
the  length  of  time  when  one  service  station  is  operating  and  the  length 
of  time  when  both  service  stations  are  operating  simultaneously  during 
the  busy  period,  respectively.  Then,  we  have 


B - B^  + B2.  (3.5) 

Hence,  from  relation  (3.4)  and  (3.5),  Bj^  and  B2  can  also  be  represented 
in  terms  of  Ug,  U2 , and  for  all  d. 


Lemma  3 . 1 

i)  The  length  of  time  when  exactly  one  service  station  is  operating 
during  the  busy  period,  Bj^  can  be  expressed  by 


Ur 


+ U2 


D-1 

E 

d-1 


Ui  + E_  Y^ 


+ U. 


for  D = 0 
for  D - 1 

for  D > 2. 


(3.6) 


ii)  The  length  of  time  when  both  service  stations  are  operating 
simultaneously  during  the  busy  period  can  be  expressed  by 
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V 


D 

d«l 


for  D = 0 
for  D > 1 


(3.7) 


where  V is  the  random  variable  with 

P[V=0]  = 1. 

Next,  we  will  characterize  Uq,  U^,  U2 , and  based  on  the 
fluctuation  of  the  number  of  customers  in  the  system  in  order  to 
characterize  the  busy  period. 


3.4  Characterizations  of  U^,  IJ.,  , , Y.|  and  V.| 

During  a busy  period,  the  number  of  customers  in  the  system  is 
increased  by  1 or  decreased  by  1 instantly  due  to  a transition  which 
could  be  a new  customer's  arrival  or  a customer's  departure  from  the 
system.  However,  after  the  number  of  customers  in  the  system  reaches  x 
during  the  busy  period,  the  system  state  x sojourns  for  a random  length 
of  time  until  it  becomes  (x+1)  or  (x-1) . Thus,  the  length  of  time  from 
the  instant  the  number  of  customers  in  the  system  reaches  x to  the  first 
time  the  number  of  customers  in  the  system  reaches  y can  be 
characterized  by  both  (i)  the  random  length  of  sojourn  time  between  two 
successive  transitions  and  (ii)  the  random  number  of  transitions  until 
the  number  of  customers  in  the  system  reaches  y for  the  first  time.  To 
do  so,  define  the  following  random  variables: 

1)  U ; the  length  of  time  from  the  instant  the  number  of  customers  in 
the  system  reaches  x to  the  first  time  the  number  of  customers 
in  the  system  reaches  y during  a busy  period. 
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2)  L ; the  number  of  transitions  during  U where  L is  a nonnegative 

integer, 

3)  : the  length  of  sojourn  time  between  the  (n-l)th  and  the  nth 

transitions  during  U where  n-l,2,---,L, 

4)  J (n)  ; the  number  of  customers  in  the  system  just  after  the  nth 

transition  during  U where  n-1,2, • • • ,L. 


J (n) 


Figure  3.2  The  fluctuation  of  the  number  of  customers 
in  the  system  during  U 

Hence,  U can  be  expressed  in  terms  of  ( n-l,2,-'-,L  ) such  as 

U - Tf  + T2  + T3  + •••  + + T^.  (3.8) 

In  a Markov  process,  every  time  epoch  at  which  a transition  takes 
place  is  a regeneration  point  which  implies  that  development  of  future 
process  is  independent  of  the  history  of  process  (see  Prabhu[20]). 
Hence,  the  T^'s  are  identically  distributed  and  independent  of  L. 
Therefore,  for  any  given  system  states  x and  y,  U can  be  characterized 
from  relation  (3.8)  if  the  random  variables  T^  and  L are  evaluated. 
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3.4.1  Evaluation  of  the  Length  of  Time 

Between  Two  Successive  Transitions  during  U 

Using  the  following  theorem,  the  random  variable  for  all  n can 
be  evaluated. 

Theorem  3 . 3 

1)  If  one  service  station  is  operating,  then  is  exponentially 
distributed  with  parameter  m(1+p)- 

2)  If  both  service  stations  are  operating  simultaneously,  then  is 
exponentially  distributed  with  parameter  /i(2+p). 

Proof  : From  the  basic  assumptions  of  system  operation  in  Chapter  1, 

the  random  variables  and  $2  are  exponentially  and  independently 

distributed  with  parameter  A,  n and  n,  respectively. 

1)  Since  - min  { A2^,Sj^},  we  have 

P[T^>t]  = P[  min{A2^,S]^}  > t] 

- P[A^^>t  , S^>t] 

- P[A^>t]  P[S^>t] 

- e-^t  g-^t 

_ g-(A+;i)t_ 

2)  Similarly,  using  - min  { A^,Sj^,S2  ) yields 

P[T3^>t]  - P[  min{A^,S3^,S2)  > t] 

- P[A^>t  , Sj^>t  , S2>t] 

- P[Aj^>t]  P[Sj^>t]  P[S2>t] 

- e-^*"  e-^^ 

= e‘ 


Thus,  the  theorem  follows  immediately. 


Q.E.D. 
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3.4.2  Evaluation  of  the  Number  of  Transitions  during  U 

In  the  stochastic  process  {J(n),  n-l,2,---,L),  the  fluctuations  of 
the  number  of  customers  in  the  system  is  a random  walk.  Hence,  L,  the 
number  of  transitions  during  U is  dictated  by  the  one  step  transition 
probabilities  given  the  initial  and  the  final  levels  of  the  number  of 
customers  in  the  systems.  Let  the  one  step  transition  probability 
matrix  of  J(n)  during  U be  H and  also  let 
p - P[J(n+l)-j+l|j(n)-j] 

and 

q - P[J(n+l)-j-l|j(n)-j] , 

where  state  j is  not  the  final  level  of  the  number  of  customers  in  the 
system,  y during  U. 

Since 

P + q = 1,  (3.9) 

we  have 

P[J(n+l)-j |j(n)-j]  - 0. 

If  state  j is  the  final  level  of  the  number  of  customers  in  the  system, 
y which  forms  an  absorbing  barrier,  then 

p - P[J(n+l)-y+l|j(n)-y]  - 0, 
q - P[J(n+l)-y-l|j(n)-y]  - 0, 

P[J(n+l)-y| J(n)-y]  - 1. 

In  general,  H exhibits  a form  with  an  absorbing  barrier  at  y if  we 
assume  without  loss  of  generality  x>y,  appearing  in  Figure  3.3. 
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Figure  3.3  H,  the  one  step  transition  probability 
matrix  of  the  number  of  customers 
in  the  system  during  U 
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The  values  of  p and  q are  obtained  using  the  following  theorem. 


Theorem  3.4 


1)  If  one  service  station  is  operating,  then 


and  q 


(1+p)  (1+p) 

2)  If  both  service  stations  are  operating  simultaneously,  then 


(3.10) 


and  q — 


(2+p)  (2+p) 

Proof  : 1)  By  Cinlar  [7]  and  Taylor  and  Karlin[24],  or 


(3.11) 


p - P[A3^<Sj^]  - 


Xe  dxdy 


0<x<y<a 

A 

(X+n) 

P 


a+p) 

2)  Similiarly, 

p - P[A3^<S^<S2]  + P[Aj^<S2<S^] 


Ae  '^^^e  dxdydz  + 


0<x<y<z<° 

A 


Ae  dxdzdy 


0<x<z<y<a 


(X+2fi) 

P 


2+p 


From  the  above  p's  and  relation  (3.9),  the  theorem  follows  immediately. 


Q.E.D. 
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Thus,  the  random  variable  L can  be  evaluated  by  the  one  step  transition 

probability  matrix  of  J (n) , H appearing  in  Figure  3.3  with  the  values  of 

p and  q shown  in  (3.10)  and  (3.11). 

Next,  we  will  characterize  Uq,  U2 , and  based  on  the 

results  of  evaluations  of  T and  L. 

n 

3.4.3  Characterizations  of  U^,  U.,  , and  Y-, 

Let  the  nonnegative  random  variables  Lq,  L2  and  L3  be  the 
number  of  transitions  during  Uq,  U2  and  for  all  d,  respectively. 
Then,  similar  to  U,  we  have 


'Uq  = Ti  + T2  + T3  + • ■ • + Tlq,  (3.12) 
Ui  - + T2  + T3  + •••  + Tl^,  (3.13) 
U2  = T3  + T2  + T3  + ...  + Tl2>  (3.14) 
Yi  - T3  + T2  + T3  + ...  + TL3.  (3.15) 


During  Uq,  Uj^,  U2  and  Y^,  one  service  station  is  operating.  Hence, 
the  Tj^'s  for  all  n are  identically,  independently  and  exponentially 
•distributed  with  parameter  . After  the  number  of  customers  in  the 

system  reaches  N when  both  service  stations  remain  inoperative 
temporarily,  the  number  of  customers  in  the  system  might  fluctuate 
between  1 and  M-1  before  it  reaches  0 or  M for  the  first  time.  Hence, 
one  of  the  two  possible  cases  will  occur  whether  the  number  of  customers 
in  the  system  will  reach  0 for  the  first  time  without  reaching  M which 
leads  to  Uq,  or  M for  the  first  time  without  reaching  0 occurs  which 
leads  to  Uj^  after  the  number  of  customers  in  the  system  reaches  N at  the 
beginning  of  the  busy  period.  The  initial  levels  of  the  number  of 
customers  in  the  system  in  Uq  and  U^  are  the  same,  namely  N and  the 
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final  levels  in  Uq  and  are  0 and  M,  respectively.  Thus,  Lq  and  , 
the  number  of  transitions  during  Uq  and  are  dictated  by  the  one  step 
transition  probability  matrix  of  the  niomber  of  customers  in  the  system 
as  shown  in  Figure  3.4. 

Similarly,  the  one  step  transition  probability  matrix  of  J(n) 
during  U2  and  which  governs  L2  and  , the  number  of  transitions 
during  U2  and  can  be  identified.  After  the  number  of  customers  in 
the  system  reaches  K for  the  first  time  while  both  service  stations  are 
operating  simultaneously,  the  number  of  customers  in  the  system  will 
reach  0 for  the  first  time  without  reaching  M which  leads  to  U2,  or  M 
for  the  first  time  without  reaching  0 which  leads  to  Yj^.  The  initial 
levels  of  the  number  of  customers  in  the  system  in  U2  and  Yj^  are  the 
same,  namely  K.  But  the  final  levels  in  U2  and  Y^  are  0 and  M, 
respectively.  Thus,  L2  and  , the  number  of  transitions  during  U2  and 
Y]^  are  also  dictated  by  the  same  one  step  transition  probability  matrix 
of  J(n)  during  Uq  and  U^^  which  governs  Lq  and  in  Figure  3.4. 

Comparing  Uq  and  U2,  both  cases  have  the  same  final  level  of  the 
number  of  customers  in  the  system,  namely  0 which  is  an  absorbing 
barrier.  Lq  and  L2 , the  number  of  transitions  during  U2  and  Y^^  are 
dictated  by  the  same  one  step  transition  probability  matrix  of  J(n) 
appearing  in  Figure  3.4.  The  only  difference  is  the  initial  levels  of 
the  number  of  customers  in  the  system,  i.e.,  N in  Uq  and  K in  U2 . This 
relation  between  Uq  and  U2  could  be  helpful  in  finding  certain 
characteristics  related  to  Uq  and  U2  from  each  other.  Note  that  the 
same  relation  between  U^^  and  Yj^  also  holds. 
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Figure  3.4  The  one  step  transition  probability  matrix 
of  the  number  of  customers  in  the  system 
during  Uq  and  where  the  values  of  p and  q 
are  shown  in  (3.10) 
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3.4.4  Characterization  of  V., 

Define  the  nonnegative  integer  random  variable  as  the  number  of 
transitions  during  Then,  can  also  be  represented  in  terms  of 

and  as 

Vi  - + l2  + T3  + ...  + Tl^.  (3.16) 

During  both  service  stations  are  operating  simultaneously. 
Hence,  the  T^'s  in  relation  (3.16)  are  distributed  identically, 
independently  and  exponentially  with  paramenter  /i(2+p) . After  the 
number  of  customers  in  the  system  reaches  M while  one  service  station  is 
operating,  the  number  of  customers  in  the  system  might  fluctuate  between 
K+1  and  an  indefinite  level  which  is  greater  than  or  equal  to  M.  Then, 
eventually  the  number  of  customers  in  the  system  reaches  K for  the  first 
time  after  a finite  number  of  transitions  because  the  system  is 
operating  under  steady- state  conditions.  Thus,  the  one  step  transition 
probability  matrix  of  J (n)  during  which  governs  is  given  by  Figure 
3.5. 

Based  on  the  results  of  characterizing  Uq,  , U2 , and  it  is 
possible  to  use  relation  (3.4)  to  proceed  towards  the  objectives  of  this 
research. 

In  the  next  chapter,  we  will  derive  the  probability  mass  function 
of  the  number  of  activations  or  removals  of  service  stations  during  the 
busy  using  the  results  of  this  chapter. 
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Figure  3.5  The  one  step  transition  probability  matrix 

of  the  number  of  customers  in  the  system  during 
where  the  values  of  p and  q are  shown  in  (3.11) 


CHAPTER  4 


PROBABILITY  MASS  FUNCTION  OF  THE  NUMBER  OF  ACTIVATIONS 
OR  REMOVALS  OF  SERVICE  STATIONS  DURING  THE  BUSY  PERIOD 

4 . 1 Introduction 

The  main  objective  of  this  chapter  is  to  derive  the  probability 
mass  function  of  the  number  of  activations  F or  removals  F of  service 

di  IT 

stations  during  the  busy  period.  The  explicit  form  of  the  probability 
mass  function  of  F^  or  F^  provides  opportunities  (1)  to  obtain  the 
Laplace  transforms  of  the  probability  density  functions  of  the  busy 
period,  of  the  busy  cycle,  of  the  length  of  time  when  one  service 
station  is  operating  and  of  the  length  of  time  when  both  service 
stations  are  operating  simultaneously  during  a busy  period  and  (2)  to 
develop  alternative  methodologies  i)  to  find  some  of  the  system 
characteristics  derived  in  Chapter  2 without  using  P(i,j)  and  ii) 
important  system  characteristics  which  may  not  be  obtainable  from 

In  order  to  find  the  probability  mass  function  of  F or  F 

cL  1C 

explicitly,  we  use  the  results  of  Chapter  3 and  the  properties  of  the 
classical  gambler's  ruin  problem. 
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^ • 2 Derivation  of  the  Probability  Mass  Function  of  Activation  or 


In  Chapter 

3 , we 

showed  that 

F^  = D + 1 

and  Fj.  = 

D + 

1. 

Hence,  relation 

(3.4) 

can  be  rewritten 

as 

Uq 

for 

1 — 1 

B - 

Ui  + 

Vi  -H  U2 

for 

Fa  = 2 

^a-^  f’a-2 

Ui  + 

S Vj,  + S 
r-1  r-1 

Yr 

-1-  U2 

for  F^ 

(4.1) 


(4.2) 


From  relation  (4.2),  it  is  evident  that  there  is  a one-to-one 
correspondence  between  the  value  of  and  the  distribution  of  the  busy 
period,  B.  In  other  words,  for  any  given  value  of  F , the  distribution 
of  B is  uniquely  determined  in  terms  of  Uq,  U^,  U2 , Y^,  and  for  all  r. 
This  implies  that  the  probability  that  the  number  of  activations  or 
removals  of  service  stations  during  the  busy  period  is  r is  equal  to  the 
probability  that  the  uniquely  determined  distribution  of  the  busy  period 
given  that  F^  is  r will  occur.  Hence,  we  define  the  following 
probabilities . 

1)  ” r],  r-1,2,---  ; the  probability  that  the  number  of 
activations  of  service  stations  during  the  busy  period  is  r,  i.e. , 
the  probability  mass  function  of  F or  F . 

2)  P(Uq)  : the  probability  that  the  number  of  customers  in  the  system 

will  reach  0 for  the  first  time  without  reaching  M after  the  number 
of  customers  in  the  system  reaches  N at  the  beginning  of  the  busy 
period,  or  the  probability  that  the  number  of  customers  in  the 
system  is  always  less  than  M during  the  entire  busy  period. 
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Equivalently,  this  is  the  probability  that  Uq  will  occur  after  the 
number  of  customers  in  the  system  reaches  N at  the  beginning  of  the 
busy  period. 

3)  ; the  probability  that  the  number  of  customers  in  the  system 
will  reach  M for  the  first  time  without  reaching  0 after  the  number 
of  customers  in  the  system  reaches  N at  the  beginning  of  the  busy 
period.  Equivalently,  this  is  the  probability  that  will  occur 
after  the  number  of  customers  in  the  system  reaches  N at  the 
beginning  of  the  busy  period. 

4)  P(U2)  ; the  probability  that  the  number  of  customers  in  the  system 
will  reach  0 for  the  first  time  without  reaching  M after  the  number 
of  customers  in  the  system  reaches  K when  both  service  stations  are 
operating  simultaneously  during  the  busy  period.  Equivalently, 
this  is  the  probability  that  U2  will  occur  after  the  number  of 
customers  in  the  system  reaches  K when  both  service  station  is 
operating  simultaneously  during  the  busy  period. 

5)  ; the  probability  that  the  number  of  customers  in  the  system 
will  reach  M for  the  first  time  without  reaching  0 after  the  nximber 
of  customers  in  the  system  reaches  K when  both  service  stations  are 
operating  during  the  busy  period.  Equivalently,  this  is  the 
probability  that  Y^.  for  all  r will  occur  after  the  number  of 
customers  in  the  system  reaches  K when  both  service  station  is 
operating  simultaneously  during  the  busy  period. 

6)  PCVj^)  : the  probability  that  the  number  of  customers  in  the  system 

will  reach  eventually  K for  the  first  time  after  the  number  of 
customers  in  the  system  reaches  M when  one  service  station  is 
operating  during  the  busy  period. 
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Hence,  it  is  easy  to  deduce  the  following  relations  from  the  definitions 
of  Uq , U^,  U2 , and  for  all  d under  steady- state  conditions; 


P(Uq) 

+ P(U^)  - 1, 

(4.3) 

P(U2) 

+ P(Yj^)  - 1, 

(4.4) 

P(vp 

- 1. 

(4.5) 

Furthermore,  according  to  relation  (4.2),  the  following  theorem  may  be 

stated  based  on  the  one-to-one  correspondence  between  the  value  of  F 

and  the  uniquely  determined  distribution  of  the  busy  period  given  the 

value  of  F„ . 

a 

Theorem  4.1 

The  probability  mass  function  of  F.  or  F^,  P[F,-r]  for  r-1,2,3,--- 

Si  iL  Si 

can  be  obtained  from  the  following  relation: 


P[F^-r] 


P(Uq)  for  r - 1 


P(U3^)P(U2)[P(V3^)]'^-1[P(Y^)]’^-2  for  r > 2. 


(4.6) 


or  using  relation  (4.4),  relation  (4.6)  is  equivalent  to 


P[Fa-r] 


P(Uq)  for  r - 1 


(4.7) 


P(U^^)P(U2)  [P(Y^)] 


r-2 


for  r > 2. 


Theorem  4.1  provides  the  fundamental  relation  to  obtain  PfF  -rl  for 
all  r;  this  depends  only  on  P(Uq) , P(U^) , P(U2)  and  P(Y^) . 

Our  next  objective  is  to  derive  explicit  forms  for  the  above  four 
probabilities  in  terms  of  A,  fx,  K,  N and  M in  order  to  obtain  P[F  =r] 
for  all  r explicitly.  To  do  so,  we  use  the  properties  of  the  gambler's 
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ruin  problem  because  of  the  similarity  between  the  fluctuations  of  the 
number  of  customers  in  the  system  in  a Markovian  death  & birth  process 
and  the  duration  of  a game  in  the  gambler's  ruin  problem. 


4.2.1  Evaluations  of  P(U(^)  and  P(U-|) 

Once  the  busy  period  is  initiated  when  the  number  of  customers  in 
the  system  reaches  N while  both  service  stations  remain  inoperative 
temporarily,  the  number  of  customers  in  the  system  will  reach  0 for  the 
first  time  without  reaching  M with  probability  P(Uq)  or,  M for  the  first 
time  without  reaching  0 with  probability  P(U^) . However,  the 
fluctuation  of  the  number  of  customers  in  the  system  during  Uq  and  is 
dictated  by  the  one  step  transition  probability  matrix  of  the  nxamber  of 
customers  in  the  system  as  exhibited  in  Figure  3.4.  During  Uq  and  U^, 
the  number  of  operating  service  stations  is  one.  The  values  of  p and  q, 
the  elements  of  the  one  step  transition  probability  matrix  of  J(n)  are 
given  by  (3.10)  such  as 

P 1 

P and  q . (4.8) 

1+P  1+p 

In  order  to  evaluate  P(Uq) , let 

$q(N)  - P(Uq).  (4.9) 

Then,  $q(N)  satisfies  the  following  second  order  difference  equations: 


$q(N)  = p $q(N+1)  + q 'I'q(N-I)  1 < N < M-1 

with  boundary  conditions 


$o(0) 

$q(M) 


(4.10) 


(4.11) 


0. 


(4.12) 
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The  solution  of  the  difference  equation  above 


(q/p)^  - (q/p)^ 
(q/p)^  - 1 


for  p 7^  q 


M-N 


M 


for  p - q. 


(4.13) 


Substituting  the  values  of  p and  q in  (3.10)  or  (4.8)  and  simplifying 
yield 


P(Uq)  - $o<N)  - 


1-p 


M-N 


1-p 

M-N 

M 


M 


for  p / 1 


for  p ” 1. 


(4.14) 


Note  that  p q and  p =•  q are  equivalent  to  p ^ 1 and  p = 1. 

P(U]^)  can  be  obtained  from  either  solving  a difference  equation 
similar  to  (4.10)  or  using  relation  (4.3)  and  the  results  of  P(Uq)  in 
(4.14).  Hence,  we  obtain 


P(up 


M-N  M 
P -P 


1-p 


M 


N 

M 


for  p ^ 1 


for  p ” 1. 


4.2.2  Evaluations  of  P(U^)  and  P(Y-^ ) 


(4.15) 


Suppose  the  number  of  customers  in  the  system  reaches  K when  both 
service  stations  are  operating  simultaneously.  Then,  the  number  of 
customers  in  the  system  will  reach  0 for  the  first  time  without  reaching 
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M with  probability  P(U2) , or  M for  the  first  time  without  reaching  0 
with  probability  P(Yj^)  . We  may  construct  and  solve  in  a similar  way  a 
second  order  difference  equation  to  obtain  P(U2>  and  P(Yj^)  . However,  we 
can  deduce  P(U2>  and  P(Yj^)  directly  from  the  results  of  P(Uq)  and  P(U^) 
in  (4.14)  and  (4.15).  One  service  station  is  operating  during  Uq  and  U2 
and  both  Uq  and  and  U2  are  first  passage  times  to  reach  0 without 
reaching  M.  This  leads  to  the  same  values  of  p and  q as  in  (4.8),  which 
are  the  elements  of  the  one  step  transition  probabilities  of  J(n).  The 
only  difference  between  Uq  and  U2  is  the  initial  levels  of  the  number  of 
customers  in  the  system,  which  are  N in  Uq  and  K in  U2 . Hence,  P(U2) 
are  obtained  directly  from  the  results  of  P(Uq)  by  replacing  N by  K. 
Similar  relation  also  holds  between  Uj^  and  Y^^.  Thus,  using  the  results 
of  P(Uq)  and  P(U^)  in  (4.14)  and  (4.15),  we  replace  N by  K to  obtain 


P(U2) 


1-p 


M-K 


1-p 

M-K 

M 


M 


for  p / I 


for  p - 1, 


(4.16) 


P(Yi) 


.M-K 

p -p 


1-p 


M 


K 

M 


for  p / 1 


for  p - 1. 


(4.17) 


Note  that  because  of  the  similarities  between  the  Markovian  birth 
and  death  process  and  the  classical  gambler's  ruin  problem,  P(Uq),  P(U2) 
P(U2^  and  P(Y^)  can  also  be  deduced  from  the  gambler's  ultimate  ruin 
probabilities  (see  Feller[10],  pp  344). 
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Since  all  four  probabilities  P(Uq)  , P(U^)  , P(U2>  and  P(Yj^)  are 
derived,  the  probability  mass  function  of  the  number  of  activations  or 
removals  during  a busy  period,  P[F^-r]  for  r-1,2,3,---  can  be  obtained 
using  relation  (4.7)  and  the  explicit  values  of  P(Uq) , P(U^) , P(U2>  and 
P(Y3^). 

4.2.3  Derivation  of  the  Probability  Mass  Function  of 

Substituting  the  explicit  results  of  P(Uq)  , P(Uj^) , P(U2>  and  P(Yj^) 
in  (4.15)  to  (4.18)  into  relation  (4.7)  yields  the  probability  mass 
function  of  F^ , P[F,-r]  for  all  r as 

cL  a 


for  p ^ 1 


P[F^-r] 


1-p 


M-N 


1-p 


M 


for  r - 1 


(l-p^-'<)(p‘^-N.,M) 

(1-P^)2 


M-K  M 

p -p 


r-2 


1-p 


M 


(4.18) 


for  r > 2, 


for  p - 1 


P[F^-r] 


M-N 

M 

N(M-K) 


for  r = 1 


K 


M 


r-2 


(4.19) 


for  r > 2. 
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From  relation  (3.2)  and  P[F^-r]  in  (4.18)  and  (4.19),  the 
probability  mass  function  of  D,  denoted  by  P[D  = d]  for  d=  0,1,2,- ••  is 
also  obtained  as 


for  p 7^  1 


P[D-d] 


1-p 


M-N 


1-p 


M 


for  d - 0 


(l-pM-K)(pM-N.pM) 


P -P 


d-1 


(4.20) 


1-p 


M 


for  d > 1, 


for  p - 1 


P[D-d] 


M-N 

M 

N(M-K) 


for  d - 0 


K 

M 


d-1 


for  d > 1. 


It  is  easy  to  verify  that 


(4.21) 


CO  00 

S P[F^-r]  - S P[D-d]  - 1.  (4.22) 

r-1  d-0 

Hence,  P[F^— r]  and  P[D—d]  in  (4.18)  to  (4.21)  are  proper  probability 

mass  functions  of  F and  D. 

a 

In  the  next  chapter,  we  will  derive  the  Laplace  transforms  of  the 
probability  mass  function  of  the  busy  period  using  the  result  of  the 
probability  mass  function  of  the  niamber  of  activations  or  removals  of 
service  stations  during  the  busy  period. 


CHAPTER  5 


LAPLACE  TRANSFORM  OF  THE  PROBABILITY  DENSITY  FUNCTION 
OF  THE  BUSY  PERIOD 

5 . 1 Introduction 

According  to  the  results  characterizing  the  random  variables  Uq, 

U2 , Y^  and  obtained  in  Chapter  3,  each  random  variable  depends  on 
the  number  of  transitions  and  the  length  of  time  between  two  successive 
transitions.  The  length  of  times  between  two  successive  transitions  are 
identically,  independently  and  exponentially  distributed  and  independent 
of  the  number  of  transitions.  Thus,  it  is  possible  to  obtain  the 
probability  density  functions  of  Uq,  U2 , Yj^  and  if  the 

probability  mass  functions  of  the  the  number  of  transitions  during  Uq , 
Uf,  U2,  Y^  and  respectively,  are  derived. 

In  this  chapter,  we  first  derive  the  probability  generating 
functions  of  the  number  of  transitions  during  Uq,  Uj^,  U2 , Y^  and 
Then,  using  these  results,  we  proceed  towards  finding  the  Laplace 
transforms  of  the  probability  density  functions  of  Uq,  Uj^  , U2 , Y^^  and 
which  are  the  basic  components  of  the  length  of  times  B,  C,  B^^  and  B2 
related  to  the  busy  period.  Finally,  we  derive  the  Laplace  transforms 
of  the  probability  density  functions  of  B,  C,  B^^  and  B2.  The  results  of 
this  chapter  will  be  used  in  Chapter  6 to  find  certain  system 
characteristics  without  using  the  joint  probability  mass  function 
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5 . 2 Probability  Generating  Functions  of  the  Number  of  Transitions 

during  U,^,  U-,  , U^,  Y-,  and 

In  order  to  find  the  probability  generating  functions  of  the  number 
of  transitions  during  Uq,  U2,  and  we  also  use  the  properties 
of  the  classical  gambler's  ruin  problem.  Hence,  let  P[Lq  - m]  for 
n>”0,l,2,---  be  the  probability  that  the  number  of  transitions  during  Uq 
is  m.  Also  let  Gq(N,z),  Gj^(N,z),  Gq(K,z),  G]^(K,z)  and  G2(M,z)  be  the 
probability  generating  functions  of  Lq,  L2 , L3  and  L^,  the  number  of 

transitions  during  Uq,  U^,  U2,  and  respectively.  For  example, 

CO 

Gq(N,z)  - 2 z“P[LQ-m].  (5.1) 

m-0 

In  order  to  find  the  probability  generating  functions  of  the  1^'s, 
we  extensively  use  the  probability  generating  function  technique  for  the 
duration  of  a game,  developed  by  Feller [10]. 

5.2.1  Gq(N,z)  and  G-,  (N,z) 

Suppose  that  the  number  of  customers  in  the  system  reaches  N when 
both  service  stations  remain  inoperative  temporarily.  Then,  the  number 
of  customers  in  the  system  fluctuates  between  1 and  M-1  before  it  will 
reach  0 for  the  first  time  without  reaching  M with  probability  P(Uq),  or 
M for  the  first  time  without  reaching  0 with  probability  P(U^) . The 
explicit  values  of  P(Uq)  and  P(Uj^)  were  already  computed  in  (4.14)  and 
(4.15)  and  given  by 


1./-N 

P(Uq) 

1-p^  ’ 

M-N  M 
P -P 

P(Ui) 

I-P^ 

(5.3) 
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Note  that 

P(Uq)  + P(U^)  = 1.  (5.4) 

In  order  to  compute  Gq(N,z),  let  PQ(N,m)  for  be  the 

unconditional  probability  that  the  number  of  transitions  is  m from  the 
epoch  when  the  number  of  customers  in  the  system  reaches  N at  the 
beginning  of  the  busy  period  to  the  epoch  when  this  number  reaches  0 
without  reaching  M.  Then,  we  have 

00 

S Po(N.m)  - PCUq),  (5.5) 

m-0 

PgCN.m) 

P[Lq  = m]  = for  all  m.  (5.6) 

P(Uq) 

The  fluctuation  of  the  number  of  customers  in  the  system  is  a 
random  walk.  Hence,  in  order  to  find  Gq(N,z),  we  use  the  methodology 
developed  by  Feller[10]  for  finding  the  probability  generating  function 
of  the  duration  of  a game  in  a classical  gambler's  ruin  problem. 

After  the  first  transition  at  the  beginning  of  a busy  period,  the 
number  of  customers  in  the  system  is  N+1  with  probability  p,  or  N-1  with 
probability  q.  Thus,  PQ(N,m)  satisfies  the  following  second  order 
difference  equation; 


PQ(N,m+l)  - p PQ(N+l,m)  + q PQ(N-l,m)  for  m > 0 and  0 < N < M (5.7) 
with  boundary  conditions 

Po(0,m)  = PQ(M,m)  - 0 for  m > 1 (5.8) 

Pq(N,0)  = 0 for  0 < N < M (5.9) 

Po(0,0)  = 1 (5.10) 

where  the  values  of  p and  q are  shown  in  (3.10). 
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Note  that  the  fluctuation  of  the  number  of  customers  in  the  system 
when  one  service  station  is  operating  during  the  busy  period  is  dictated 
by  the  one  step  transition  probability  matrix  of  the  number  of  customers 
in  the  system  appearing  in  Figure  3.4. 

Define 

CO 

go(N,z)  = 2 z™  PgCN.m).  (5.11) 

m=0 

Multiplying  (5.7)  by  z*"^^  and  adding  over  all  possible  values  of  m,  we 
have 

go(N,z)  - pz  gQ(N+l,z)  + qz  gQ(N-l,z)  for  0 < N < M (5.12) 
with  boundary  conditions 

gQ(0,z)  = 1 (5.13) 

go(M,z)  = 0.  (5.14) 


Solving  the  second  order  difference  equation  (5.12)  with  boundary 
conditions  (5.13)  and  (5.14)  and  substituting  the  values  of  p and  q in 
(3.10),  we  obtain  gQ(N,z)  as 


where 


go(N.z) 


(2z)N{  ^^(z)^-N  - ) 

{ - (?2(z)^  ) 


= 1+p  + { (l+p)^  - 4pz^ 

^2(2)  = 1+p  - { (1+p)^  - 4pz^  . 


(5.15) 


(5.16) 

(5.17) 


However,  from  relation  (5.1)  and  (5.5),  we  have 

go(N,z) 

p(Uq)  ■ 


Gq(N,z) 


(5.18) 


Thus,  substituting  gQ(N,z)  in  (5.15)  into  relation  (5.18),  we  obtain 
Gq(N,z),  the  probability  generating  functions  of  the  number  of 
transitions  during  Uq  as 


Gq(N,z)  = 


(1-/'N){  ^i(z)^  - 82(z)^  ) 


(5.19) 


Similarly,  we  can  obtain  Gj^(N,z),  the  probability  generating 
function  of  the  number  of  transitions  during  as 


5.2.2  Gq(K,z)  and  G-](K,z) 

Suppose  that  the  number  of  customers  in  the  system  reaches  K when 
both  service  stations  are  operating  simultaneously.  The  number  of 
customers  in  the  system  will  fluctuate  until  it  reaches  0 for  the  first 
time  without  reaching  M with  probability  P(U2) , or  M for  the  first  time 
without  reaching  0 with  probability  P(Y2^) . However,  the  fluctuations  of 
the  number  of  customers  in  the  system  during  U2  and  are  also  dictated 
by  the  one  step  transition  probability  matrix  of  the  number  of  customers 
in  the  system  as  shown  in  Figure  3.4.  The  values  of  p and  q which  are 
the  elements  of  the  one  step  transition  probability  matrix  of  J(n)  are 
given  in  (3.10)  because  one  service  station  is  operating  during  U2  and 
Y^ . Hence,  Gq(K,z)  and  Gj^(K,z)  can  be  obtained  using  similar  procedures 
when  finding  Gq(N,z). 

Gq(K,z)  and  G^(K,z)  can  also  be  deduced  directly  from  the  results 
of  Gq(N,z)  and  G^(N,z)  based  on  the  same  approaches  used  to  find  P(U2) 
and  P(Y^)  from  P(Uq)  and  P(U^)  in  Chapter  4.  Thus,  replacing  N by  K in 


Gi(N,z)  - 


(l-/)(2pz)N{  ^^(z)N  - ^2(z)”  ) 


(5.20) 


73 


the  results  of  Gq(N,z)  and  G]^(N,z)  in  (5.19)  and  (5.20),  we  obtain 
OqCK.z)  and  G]^(K,z),  the  probability  generating  functions  of  the  number 
of  transitions  during  U2  and  as 


Gq(K,z) 


G^(K,z) 


(l-p^)(2z)'^{  (?i(z)'^-K  - (?2(z) 


'2(2)^  ) 


(l-p^)(2pz)‘^{  (?i(z)*^  - ^2(2)*^  ) 


^K, 


(p^'^-p^){  - ^2(2) 


M 


where  ^^^(z)  and  ^2^^)  also  in  (5.16)  and  (5.17). 


(5.21) 


(5.22) 


5.2.3  G,(M,z) 


G2(M,z),  the  probability  generating  function  of  L^,  the  number  of 
transitions  during  V^,,  can  also  be  obtained  by  a similar  method  used  to 
find  Gq(N,z).  Let  P[L^  = m]  for  m=0,l,2,---  be  the  probability  mass 
function  of  L^.  Then,  G2(M,z)  is  defined  by 

00 

G2(M,z)  = S z™P[L^=m].  (5.23) 

m=0 

After  the  number  of  customers  in  the  system  reaches  M when  one 
service  station  is  operating,  the  level  of  customers  in  the  system  may 
fluctuate  between  K+1  and  an  indefinite  level  which  is  equal  to  or 
greater  than  M before  it  is  decreased  to  K for  the  first  time.  Because 
the  system  operates  under  steady-state  conditions,  the  probability  that 
eventually  the  number  of  customers  in  the  system  reaches  K during 
after  a finite  number  of  transitions  is  one.  However,  the  fluctuation 
of  the  number  of  customers  in  the  system  during  V^.  is  dictated  by  the 
one  step  transition  probability  matrix  of  J(n)  appearing  in  Figure  3.5. 
The  values  of  p and  q,  which  are  the  elements  of  the  one  step  transition 
probability  matrix  of  J(n)  are  shown  in  (3.11)  because  both  service 
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stations  are  operating  simultaneously  during  V^,.  Hence,  if  we  define 


P2(M,m)  = P[L^  = ra]  for  all  m,  (5.24) 

then,  we  have  the  following  second  order  difference  equation: 

P2(M,m+l)  = p P2(M+l,m)  + q P2(M-l,m)  for  m > 0 and  K+1  < M (5.25) 
with  boundary  conditions 

P2(K,m)  - 0 for  m > 1 (5.26) 

P2(K,0)  - 1 (5.27) 

P2(M,m)  - 0 for  m > 1.  (5.28) 


Multiply  (5.25)  by  and  add  over  all  possible  values  of  m.  We 

obtain 


G2(M,z)  = pz  G2(M+1,z)  + qz  G2(M-l,z)  for  K < M (5.29) 

with  boundary  condition 

G2(K,z)  - 1.  (5.30) 

In  order  to  solve  completely  the  second  order  difference  equation 
in  (5.29)  , we  need  to  specify  another  condition.  The  following  indirect 
approach  may  be  used  to  determine  G2(M,z). 

As  shown  in  Figure  3.5,  the  one  step  transition  probability  matrix 
of  the  fluctuation  of  the  number  of  customers  in  the  system  during  V^, 
has  an  infinite  number  of  states  because  the  number  of  customers  in  the 
system  can  be  increased  to  infinity.  This  situation  eliminates  one 
boundary  condition  of  the  difference  equation  (5.29)  . Hence,  consider  a 
problem  in  which  we  impose  an  upper  restriction  in  the  fluctuation  of 
the  number  of  customers  in  the  system,  say  X where  X < M.  Then,  the  one 
step  transition  probability  matrix  for  the  fluctuation  of  the  number  of 
customers  in  the  system  will  have  a finite  number  of  states.  Thus,  the 
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number  of  customers  in  the  system  can  only  fluctuate  between  K+1  and  X-1 
before  it  reaches  K for  the  first  time  without  reaching  X,  or  X for  the 
first  time  without  reaching  K.  This  implies  that  we  can  have  two 
boundary  conditions  for  the  second  order  difference  equation  in  (5.29). 

solution  of  the  following  second  order  difference 

equation: 


g2(M,z)  - pz  g2(M+l,z)  + qz  g2(M-l,z) 

for  K < M 

(5.31) 

with  boundary  conditions 

g2(M,z)  - 1 

(5.32) 

g2(X,z)  = 0. 

(5.33) 

Suppose  that  we  increase  the  value  of  X to  infinity.  Then,  the  number 
of  customers  in  the  system  can  fluctuate  without  any  upper  restrictions 
until  it  is  absorbed  to  K.  Therefore,  as  the  value  of  X goes  to 
infinity,  g2(M,z)  becomes  G2(M,z)  such  that 


G2(M,z)  = lim  g2(M,z).  (5.34) 

X->« 

Hence,  we  obtain  the  probability  generating  function  of  the  number  of 
transitions  during  G2(M,z)  as 


G2(M,z) 


2+p  - {(2+p)2  - 8pz2)V2 

2pz 


M-K 


(5.35) 


All  these  explicit  forms  of  the  probability  generating  functions, 
Gq(N,z),  G]^(N,z),  Gq(K,z),  Gq(N,z)  and  G2(M,z)  are  the  basic  information 
necessary  to  derive  the  Laplace  transforms  of  the  probability  density 
functions  of  Uq,  , U2 , and  discussed  in  the  next  section. 
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The  results  of  the  probability  generating  functions  could  be  used 
to  derive  E[Uq]  , E[U]^],  E[U2]  , E[Yj^]  and  E[Vj_].  Let  E[Lq]  , E[L^]  , 
E[L2],  E[L3]  and  E[L^]  be  the  expected  value  of  Lq  , , L2 , L3  and  L^, 

respectively.  Then,  from  relation  (3.12)  to  (3.16),  we  have 


E[Uq]  - E[Ti]  E[Lo], 
E[U]^]  = E[T^]  E[L^]  , 
E[U2]  » E[T]_]  E[L2]  , 
E[Y^]  - E[TJ  E[L3], 
E[Vj^]  - E[T^]  E[L^]. 


E[Tj^]  was  obtained  as 


mi] 


1 

1 

p(2+p) 


for  0-1 


for  0=2. 


(5.36) 

(5.37) 

(5.38) 

(5.39) 

(5.40) 


(5.41) 

(5.42) 


When  computing  E[Uq]  , E[U]^]  , E[U2]  and  E[Y^],  E[T]^]  for  the  case  0-1 
in  (5.41)  is  considered,  while  when  computing  E[V^],  E[T]^]  for  the  case 
0 - 2 in  (5.42)  is  considered. 

E[Lq],  E[L^],  E[L2],  E[L3]  and  E[L^]  are  obtained  respectively  from 
Gq(N,z),  G]^(N,z),  Gq(K,z),  G-j^(K,z)  and  G2(M,z)  by  differentiating  with 
respect  to  z and  setting  z - 0.  We  obtain  as  an  example  E[Lq]  for 
E[Uq].  The  other  expressions  in  (5.37)  to  (5.40)  are  obtained 
similarly.  E[Lq]  is  obtained  as 


E[Ln] 


(l+p)[N(l+/-N)(l-pM)  . 2M(p^-N-/)j 
(l-p)(l-p^)(l-p^-N) 


(5.43) 


Substituting  E[Tj^]  in  (5.41)  and  the  above  E[Lq]  into  relation  (5.36) 
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yields 


E[Uq] 


(5.44) 


The  explicit  forms  of  E[Uq]  , E[U]^]  , E[U2]  , E[Yj^]  and  E[V^]  will  be  used 
to  develop  alternative  methodologies  to  find  certain  system 
characteristics  in  Chapter  7. 


5 . 3 Laplace  Transforms  of  the  Probability  Density 
Functions  of  U^,  , Y.|  and 


The  Laplace  transforms  of  the  probability  density  functions  of  Uq, 
Uf,  U2 , Y^^  and  may  also  be  used  for  finding  the  values  of  E[Uq]  , 
E[U^^]  , E[U2],  E[Yj^]  and  E[Vj^]  . More  importantly  however,  they  are  used 
to  find  the  Laplace  transforms  of  the  probability  density  functions  of 
B,  C,  and  B2  explicitly. 

Let 


1)  fy^(x)  for  0<x<«>  ; the  probability  density  function  of  Uq , 


Laplace  transforms 

of  the  probability  density  function  of  Uq,  U^^,  U2 , Y^^  and  V^, 
respectively,  i.e.,  for  example. 


(s) 


e'®^  fy^Cx)  dx. 


(5.45) 


In  order  to  find  fuQ(s),  fuj^(s), 

the  convolution  properties  for  the  sum  of  independent  random  variables 
and  the  results  of  Gq(N,z),  G^(N,z),  Gq(K,z),  Gj^(K,z)  and  G2(M,z) 
obtained  in  (5.19)  to  (5.22)  and  (5.35). 
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5.3.1 


(s) 


Define  the  following; 

i)  f.j(x)  for  0<x<®  ; the  common  probability  density  function  of  T^, 

the  length  of  time  between  two  successive  transitions  during  Uq 
(Note  that  the  T^'s  are  identically,  independently  and  exponentially 
distributed  with  parameter  , and  independent  of  Lq.), 

ii)  f.p(s)  ; the  Laplace  transform  of  f^(x)  such  as 


f^(s) 


fj(x)  dx, 

Jo 


(5.46) 


iii)  fu^Cx  I Lo=m)  for  all  m ; the  conditional  probability  density 

function  of  Uq  given  that  the  number  of  transitions  during  Uq  is 


m. 


Then,  we  have 


^i(l+p) 

= . (5.47) 

s+m(1+p) 

Also,  using  relation  (3.12)  yields 

“ fj*^“^(x)  for  all  m,  (5.48) 

where  f.p  -'(x)  denotes  the  m-fold  convolution  of  f^(x)  with  itself. 
Thus,  fy^(x)  is  obtained  from 

CD 

^ = "“I  • (5.49) 

m~0 

Substituting  relation  (5.48)  into  relation  (5.49)  and  then  taking  the 
Laplace  transform  on  both  sides,  we  obtain 


00 

fu  (s)  = Z [f^(s)]"'  P[Lq  = m]. 

m-0 


(5.50) 
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However,  from  relation  (5.1),  we  have 

“ GQ[N,f^(s)].  (5.51) 

Relation  (5.51)  implies  that  fy^(s)  can  be  obtained  by  substituting 
f^(s)  in  (5.47)  instead  for  z in  Gq(N,z)  in  (5.32).  Hence, 


of  the  probability  density  function  of  Uq 

is  obtained  as 

(l-p^)(2^i)N[ 

{jT]^(s)  )^' 

- J 

(5.52) 

(1-P^-N)[ 

{7T]^(S) 

- ] 

where  ’t^(s) 

“ s+A+/i  + { 

(s+A+/i)^ 

- 4Am 

(5.53) 

7T2(s) 

— s+A+/i  - { 

(s+A+p)^ 

- 4Am  )V2. 

(5.54) 

5.3.2  fy^(s) , fy^(s),  fy^(s)  and  fy^(s) 


Similarly,  we  can  obtain  fu2(s),  explicitly. 

During  , U2  and  Y^,  one  service  station  is  operating.  Hence, 
substituting  f^(s)  in  (5.47)  instead  for  z in  Gj^(N,z),  Gq(K,z)  and 
G]^(K,z)  derived  in  (5.20)  to  (5.22),  we  obtain  fu^(s),  ^2(3)  and 

Laplace  transforms  of  the  probability  density  function  of 
U^,  U2  and  as 


(1-P^)(2A)^-N[  {^^(s)}N  - {,t2(s))N  ] 

(p^'N  -p^)[  {7r]^(s))^  - {7r2(s))^  ] 

(l-p^)(2/i)^[  {:r^(s))f^->^  - ] 

(1-P^  ^)  [ (?T]^(S))^  - {^2(3))^  ] 

(l-p^)(2A)^’^[  {7T3^(s))‘^  - {7T2(s))^  ] 

(p^-^  -/)[  {7T3^(s))^  - {7T2(s)}^  ] 


(5.55) 


(5.56) 


(5.57) 


where  7Tj^(s)  and  ?T2(s)  are  in  (5.53)  and  (5.54). 
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However,  during  both  service  stations  are  operating 
simultaneously.  Thus,  the  length  of  times  between  two  successive 
transitions  are  identically,  independently  and  exponentially  distributed 
with  parameter  n(2+p) . The  Laplace  transform  of  the  common  probability 
density  function  of  the  length  of  times  between  two  successive 
transitions  during  becomes  f^(2+p) / (s+p(2+p) ) . Substituting  this 
value  instead  for  z in  G2(M,z)  in  (5.11),  we  obtain  fy^(s),  the  Laplace 
transform  of  the  probability  density  function  of  as 


s+X+2n  - ( (s+X+2^i)^  - 8A/1 


2A 


M-K 


(5.58) 


Next,  we  will  derive  the  Laplace  transforms  of  the  probability 
density  functions  of  B,  C,  and  B2 . 


5 ■'i  Laplace  Transforms  of  the  Probability  Density  Functions 

of  B,  C,  B-|  and  B^ 


the  previous  section,  we  obtained  the  basic  components  for 
deriving  the  Laplace  transforms  of  the  probability  density  functions  of 
the  busy  period  B,  of  the  busy  cycle  C,  of  the  length  of  time  when  one 
®®tvice  station  is  operating  and  of  the  length  of  time  when  both  service 
stations  are  operating  simultaneously  during  the  busy  period,  B]^  and  B2 . 
We  now  proceed  with  the  necessary  definitions  to  derive  these  Laplace 
transforms . 

Define  the  following: 

i)  fg(x) , fg^(x)  , fg^(x)  and  f|-j(x)  for  0 < x < « ; the  probability 
density  function  of  B,  Bj^,  B2  and  C,  respectively. 
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ii)  fg(s),  fg^(s),  f^Cs);  the  Laplace  transforms  of  fg(x), 

^respectively , 

iff)  fj(s)  ; the  Laplace  transform  of  the  probability  density  function 
of  the  idle  period, 

iv)  fg(x|F^-r),  fg^(x|F^=r)  and  ^ 

conditional  probability  density  function  of  B,  Bj^  and  B2  given  that 
Fa>  the  number  of  activations  or  removals  of  service  stations 
during  the  busy  period  is  r, 

■v)  fu^(x),  fy^(x),  fy^Cx),  fy(x)  and  fy^(x)  for  0 < x < ® ; the 

probability  density  functions  of  the  random  variables  U^,  U2 , Y^,  V 
and  respectively, 

; the  Laplace  transform  of  fy(x)  . 

5.4.1  Derivation  of  fg(s) 

We  first  derive  fg(s)  based  on  relation  (4.3)  using  the  convolution 
properties  of  the  sum  of  independent  random  variables.  The  result  of 
the  probability  mass  function  of  F^,  P[F^  - r]  obtained  in  (4.19)  and 
the  Laplace  transforms  of  the  probability  density  functions  of  Uq,  , 

^2’  ^1  ^1  ot)tained  in  (5.52)  to  (5.58)  are  required  for  proceeding 

to  our  objectives. 

From  relation  (4.3),  we  have  fg(x|F^-r)  for  all  r as 

= fu^(x)*fy^(x)*fu^(x),  (5.60) 

*(t-l)(x)*fY^*(t-2)(x)*f^^(x)  for  r > 3,  (5.61) 
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where  * denotes  convolution  operation  and  denotes  the  r-fold 

convolution  of  f(x)  with  itself. 

Thus,  we  obtain  fij(x)  is  from 


S ffiCxlF  -r)  P[F  - r] 


r-1 


(5.62) 


Substituting  fg(x|F^-r)  for  all  r in  (5.59)  to  (5.61)  into  (5.62),  and 
then  taking  the  Laplace  transform  on  both  sides,  we  obtain 


- fu^(s)P[F^-l]  + %,(s)fu  (s)fv  (s)  E [L  (s)fY  (s)]’^-2  P[F^-r]. 

(5.63) 

Substituting  P[F^-r]  in  (4.18)  into  (5.63)  and  simplifying  yield  fg(s), 
the  Laplace  transform  of  the  probability  density  function  of  the  busy 
period  as 


fR(s) 


where 


l-p 


M-N 


M 


l-p 


(l-p^-^)(pM-N  -pf^)fy^(s)fu2(s)fv^(s) 


(l-p^)[l-p^  - (p”-‘^-p^)fv^(s)fY^(s)] 


(l-p^)(2p)N[  - {^2(s))”'^ 

(1-pM-N)^  {^^(s))M  . {;t2(s))^  ] 


(5.64) 


(5.65) 


%^(s) 


(1-P^)(2A)M-N[  [n^(s))^  - {^2<s))^  ] 

(p^'N.pM)j  . {;t2(s))^  ] 


(5.66) 


(l-p^)(2p)^[  {,r^(s.))^-K  - ] 

(l-p”-^)[  {7T]^(S)}^  - {7T2(s))^  ] 


(5.67) 


fy]^(s) 


(l-p^)(2A)^’^-"^[  {7T3^(s)}'^  - {7T2(s)) 

(p^-^-p^)[  {7Tj^(s))^  - (7r2(s))^  ] 


vM-K 


K 


,K 


(5.68) 
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fv^(s) 


s+X+2fj,  - { (s+X+2/j)^  - 8A^ 


M-K 


2A 


(5.69) 


7T^(s)  - s+A+;i  + { (s+A+/i)2  - 4A/i  (5.70) 

7T2(s)  - s+A+m  - { (s+A+m)^  - 4A/1  (5.71) 


From  the  Laplace  transform  of  the  probability  density  function  of 
the  busy  period  obtained  in  (5.61)  to  (5.68),  we  can  deduce  the  Laplace 
transforms  of  the  probability  density  functions  of  the  busy  periods  for 
the  special  cases  of  the  M/M/2  system  operating  under  the  (0,K,N,M) 
policy. 

The  M/M/2  system  operating  under  the  (0,1, 1,2)  policy  is  equivalent 
to  the  ordinary  M/M/2  system.  Hence,  setting  K-N-1  and  M=2  in  relations 
(5.64)  to  (5.69)  yields  the  Laplace  transform  of  the  probability  density 
function  of  the  busy  period  for  the  ordinary  M/M/2  system  as 

5 n?-  ^5.72) 

s+A  + {(s+A+2/i)^  - 8Am)^/'^ 

thus,  recovering  Conolly's  result  (see  Conolly[8]). 

Another  special  case  of  our  model  is  the  (0,K,N,«>)  policy  which  is 
equivalent  to  the  M/M/1  system  operating  under  the  N policy,  originated 
by  Yadin  and  Naor[27] . Hence,  taking  the  limit  in  relation  in  (5.64)  to 
(5.69)  as  M goes  to  infinity  yields  the  Laplace  transform  of  the 
probability  density  function  of  the  busy  period  for  the  M/M/1  system 
operating  under  the  N policy  as 


2m 


iN 


s+X+n  + {(s+X+^i)^  - 


(5.73) 
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Since  the  ordinary  M/M/1  system  is  a special  case  of  the  M/M/1 
system  operating  under  the  N policy  (N-1) , the  Laplace  transform  of  the 
probability  density  function  of  the  busy  period  for  the  ordinary  M/M/1 
system  is  immediately  obtained  from  (5.73)  (see  Conolly[8]  or 
Kleinrock[ 15 ] ) . 

We  thus  deduce  that  the  busy  period  for  the  M/M/1  system  operating 
under  the  N policy  is  the  sum  of  N busy  periods  of  the  ordinary  M/M/1 
system. 

Next,  we  will  derive  fc(s),  fB^(s)  and  ^32(2),  the  Laplace 
transforms  of  the  probability  density  functions  of  C,  and  82- 

5.4.2  Derivation  of  fc(s) , fg^(s)  and  fg^C^) 


f(](s),  the  Laplace  transform  of  the  probability  density  function  of 
the  busy  cycle  can  be  obtained  by  applying  the  convolution  properties  to 
relation  (2.29).  From  relation  (2.31),  we  have 


fl(s) 


s+A 


(5.74) 


However,  from  relation  (2.29),  we  have 


fc(s)  = fj(s)  fg(s). 


(5.75) 


Thus,  substituting  fj(s)  in  (5.74)  and  the  result  of  fB(s)  in  (5.64)  to 
(5.71)  into  (5.75)  yields  fQ(s)  explicitly. 

Laplace  transforms  of  the  probability 
density  functions  of  the  length  of  time  when  one  service  station  is 
operating,  B^^  and  of  the  length  of  time  when  both  service  stations  are 
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operating  simultaneously,  B2  during  the  busy  period.  Using  the  same 
procedures  to  find  fg(s)  based  on  relation  (3.6)  yields  fg^(s)  as 


fB^(s)  = 


1-p 


M-N 


1-p 


M 


/I  M-Ks.  M-N  M^t 

(f'P  )(p  "P  ) ^ ^ ' ^U2  ^ ^ 


oMn  n .M  , M-K  t ^ , 
(1-p  )[l-p  - (p  -p  )fY^(s)J 


(5.76) 


However,  the  random  variable  V in  relation  (3.7)  was  defined  as 

P[V  - 0]  - 1.  (5.77) 

Hence , we  have 


f^(s)  - P[V  - 0]  - 1. 


(5.78) 


Therefore,  using  relation  (3.7),  we  can  obtain 


1 M-N  , M-K  „M-N,,i  , 

1-P  - (P  -p 


1-p^  - (p^'^-p^)fY^(s) 


(5.79) 


Thus,  we  can  obtain  fg^(s)  and  fg^(s)  explicitly  by  substituting 

(5.65)  to  (5.71)  into 

relation  (5.78)  and  (5.79). 

Even  though  the  Laplace  transforms  of  the  probability  density 
functions  of  B,  B^  and  B2  can  be  obtained,  it  is  difficult  to  invert 
them  in  general.  However,  the  distribution  function  of  B,  B^^  and  B2  can 
be  derived  explicitly  in  terms  of  convolutions  of  the  probability 
density  functions  of  Uq,  , U2 , and  (see  e.g.  Appendix  A). 

Hence,  it  is  possible  to  investigate  the  behavior  of  the  distribution  of 
B,  B^^  and  B2  for  selected  numerical  values  of  A,  p,  K,  N and  M.  When 
p-1 , K-N=2  and  M-3 , the  distribution  function  of  Bj^ , B2  and  B for 
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different  values  of  n are  shown  in  Figure  5.1,  5.2  and  5.3, 
respectively.  In  particular,  the  distribution  function  of  B2  has  a 
concentration  of  mass  with  magnitude  P[F„=1]  at  x=0.  This  is  due  to  the 

a 

fact  that  both  service  stations  cannot  be  operating  simultaneously  when 
the  number  of  activations  or  removals  of  service  stations  during  the 
busy  period  is  one. 

In  the  next  chapter,  we  will  develop  alternative  methodologies  to 
find  the  system  characteristics  using  the  results  of  Chapter  3,  4 and  5. 
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fs^(x) 


Figure  5.1  The  probability  density  function  of  B, 
when  A*  “ 10,  3,  1 and  0.25 
given  p-1,  K-N-2  and  M - 3 
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Figure  5.2  The  mixed  distribution  function  of 
when  /j  - 10,  3,  1 and  0.25 
given  p - 1,  K-N-2  and  M - 3 
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fg(x) 


Figure  5.3  The  probability  density  function  of  B 
when  /i  — 10,  3,  1 and  0.25 
given  p-1,  K-N-2  and  M - 3 


CHAPTER  6 


ALTERNATIVE  METHODOLOGIES  TO  FIND  SYSTEM  CHARACTERISTICS 

6 . 1 Introduction 

The  primary  objective  of  this  chapter  is  to  derive  alternative 
methodologies  to  find  certain  system  characteristics  that  will  be  used 
to  construct  a total  expected  cost  function  in  Chapter  7.  Some  of  these 
characteristics  were  obtained  in  Chapter  2.  However,  the  methodology 
does  not  use  the  explicit  form  of  the  steady- state  joint  probability 
mass  function  of  the  number  of  customers  and  the  number  of  operating 
service  stations  in  the  system,  P(i,j).  It  provides  an  alternative  way 
to  derive  these  system  characteristics  using  the  results  of  Chapter  3,  4 
and  5 . 

The  following  is  a list  of  the  system  characteristics  to  be 
obtained  most  of  which  being  first  moments: 

(1)  the  expected  number  of  activations  or  removals  of  service  stations 
during  the  busy  period,  denoted  by  E[F^]  or  E[F^] , 

(2)  the  conditional  expected  length  of  times  related  to  the  busy  period 
given  that  F^  is  r for  all  r such  as 

1)  the  conditional  expected  length  the  busy  period,  denoted  by 
E[B|F^=r] , 

2)  the  conditional  expected  length  of  the  busy  cycle,  denoted  by 
E[C|F^-r], 
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3)  the  conditional  expected  length  of  time  when  one  service  station 
is  operating  during  the  busy  period,  denoted  by  E[B]^|F^=r], 

4)  the  conditional  expected  length  of  time  when  both  service  stations 
are  operating  simultaneously  during  the  busy  period,  denoted  by 
E[B2|F^-r], 

(3)  the  expected  length  of  times  related  to  the  busy  period  such  as 

1)  the  expected  length  of  the  busy  period,  denoted  by  E[B], 

2)  the  expected  length  of  the  busy  cycle,  denoted  by  E[C] , 

3)  the  expected  length  of  time  when  i service  stations  are  operating 
during  the  busy  period,  denoted  by  E[B^^]  for  i-1,2, 

(4)  the  steady- state  probabilities  such  as 

1)  the  probability  that  i service  stations  are  operating  during  the 
busy  cycle,  denoted  by  P[0-i]  for  i-0,1,2, 

2)  the  steady- state  probability  that  no  service  stations  are 
operating  and  no  customers  are  in  the  system,  denoted  by  P(0,0), 

(5)  the  expected  number  of  operating  service  stations,  denoted  by  E[0] , 

(6)  the  expected  number  of  activations  or  removals  of  service  stations 
per  unit  time,  denoted  by  E[R^]  or  E[Rj_]  . 

6 . 2 Expected  Number  of  Activations  or  Removals 
of  Service  Stations  During  the  Busy  Period 

The  expected  number  of  activations  or  removals  of  service  stations 
during  the  busy  period,  E[F^]  or  E[F^]  is  obtained  using  the  probability 
mass  function  of  F^,  P[F^-r]  obtained  in  (4.19)  and  (4.20).  Hence,  we 


have 
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E[F^]  > E[F^] 


S r P[F^  = r] 
r=l 


M-N  M 
P -P 


1 + 


1 + 


l-p 

N 

M-K 


M-k 


for  p 1 


for  p - 1. 


(6.1) 


However,  from  relations  in  (3.3)  and  (6.1),  E[D]  is  also  obtained  as 


E[D]  - S d P[D  - d] 
d-0 


M-N  M 
P -P 


\-p 

N 

M-K 


M-k 


for  p ^ 1 


for  p = 


(6.2) 


6 . 3 Conditional  Expected  Length  of  Times  related  to  the  Busy  Period 
6.3.1  E[B|F^-r]  for  all  r 


Taking  expectation  on  both  sides  of  relation  (3.4)  or  (4.3),  we 
have  for  the  conditional  expected  length  of  busy  period  given  that  F is 

3. 

r 


E[B|F^-r] 


E[Uq]  for  r-1 

E[U^]  E[U2]  + (r-l)E[Vj^]  + (r-2)E[Vj^] 


(6.3) 

for  r>2. 


In  order  to  evaluate  the  values  of  E[Uq]  , E[U]^]  , E[U2]  , E[Yj^]  and 
E[Vj^],  differentiate 
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(5.52)  and  (5.55)  to  (5.57)  with  respect  to  s and  set  s to  zero,  we 
obtain  respectively 


E[Uq] 


N(l+p^-N)(l-pM)  . 2M(p^-N.pM) 
M(l-p)(l-p^)(l-pM-N) 

N(2M-N) 

for  p - 1, 

Sfjt 


for  p 7^  1 


(6.4) 


E[U^] 


M(l+p^)(l-pN)  - N(l+pN)(l-pM) 

M(l-p)(l-pN)(l-p”) 

(M+N) (M-N) 

for  p = 1, 

6p 


(6.5) 


E[U2l  = 


K(l+p^-‘^)(l-p^)  - 2M(p^-^-p^) 
M(l-p)(l-p^)(l-p^-‘^) 
K(2M-K) 

for  p = 1, 

6p 


for  p / 1 


(6.6) 


E[Y^] 


M(l+p^)(l-p^)  - K(l+p^)(l-p^) 
/i(l-p)(l-p^)(l-p^) 

(M+K)(M-K) 

for  p - 1, 

6p 


for  p 7^  1 


(6.7) 


M-K 

E[VJ  - . (6.8) 

4(2-p) 


Note  that  E[U2]  and  E[Yj^]  can  also  be  obtained  directly  from  the 
results  of  E[Uq]  and  E[U]^]  by  replacing  N by  K.  Also  note  that 
relations  (5.36)  to  (5.40)  withGo(N.z),  G]^(N,z),  Gq(K,z),  Gj^(K,z)  and 
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G2(M,z)  in  (5.19)  to  (5.22)  and  (5.35)  yields  the  same  E[Uq]  , E[U]^], 
E[U2],  E[Y2^]  and  E[V^^]  as  in  (6.4)  to  (6.8). 

Substitute  E[Uq]  , E[Uj^]  , E[U2],  E[Y^]  and  E[V]^]  in  (6.4)  to  (6.8) 
into  relation  (6.3)  and  simplify;  we  obtain 


N(l+p^-N)(l-pM)  . 2M(p^-N-p^^) 

m(1-p)(1-p^)(1-p^‘N) 


for  = 1 and  p / 1 


E[B|F^--r] 


N(2M-N) 

6m 


for  r = 1 and  p ~ 1 


(6.9) 


M(l-P) 


2K(l-p^) 


2M(p^-K-P^) 


(l-p‘^)(l-pM-’<)  (l-pM)(l-pM-K) 


N(l+pN) 

(1-pN) 


(r-1) 


m(1-p)(2-p) 


M(3-2p+p^)  K(3-2p+p^) 


(1-P^) 


(1-p^) 


for  r > 2 
and  p / 1 


2MK  - + (r-1)  (M-K)  (6+M+K) 

6p 


for  r > 2 and  p = 1 . 


6.3.2  E[C|F^^r]  for  all  r 

The  conditional  expected  length  of  the  busy  cycle  given  that  F is 
r,  E[C|F^-r]  can  be  obtained  from 

E[C|F^=r]  = E[I]  + E[B|F^=r],  (6.10) 

where  E[I]  is  the  expected  length  of  the  idle  period  and  is  given  by 

E[I]  = N/A. 
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6.3.3  E[Bt  |F^-r]  for  all  r 

Using  a similar  procedure  for  finding  E[BlF.,=r],  the  conditional 
expected  length  of  time  when  exactly  one  service  station  is  operating 
during  the  busy  period  given  that  is  r,  E[B]^|F^-r]  can  be  obtained. 
From  relation  (3.6),  we  have 


Thus,  substituting  E[Uq]  , E[U^]  , E[U2]  and  E[Yj^]  in  (6.4)  to  (6.7)  into 
relation  (6.11)  and  simplifying  yield 


E[Uq] 


for  r - 1 


E[BjF^-r]  - 


(6.11) 


E[Uj^]  + E[U2]  + (r-2)E[V]^]  for  r > 2. 


N(l+p^-N)(i-pM)  . 2M(p^- 


m(1-p)(1-p^)(1-p^-n 


for  r - 1 and  p / 1 


N(2M-N) 


for  r = 1 and  p = 1 


E[Bi|F^-r]  - 


(6.12) 


M(l-p)  (l-p‘^)(l-p^-^)  (l-p^)(l-pM-J^)  (1-pN) 


(r-1)  M(l+p^)  K(l+p^) 


+ 


M(l-p)  (1-p^)  (l-p*^) 


for  r > 2 and  p / 1 


2MK  - + (r-l)(M-K)(M+K) 


for  r > 2 and  p - 1. 
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6.3.4  E[B^|F^=r]  for  all  r 


From  relation  (3.7),  we  obtain  the  conditional  expected  length  of 
time  when  both  service  stations  are  operating  simultaneously  during  the 
busy  period  given  that  F^  is  r,  E[B2|F^-r]  as 


E[B2|F^-r] 


0 for  r - 1 

(r-l)E[V^]  for  r > 2. 


(6.13) 


Substituting  E[Vj^]  in  (6.8)  into  the  above  relation  yields 


E[B2|F^-r] 


0 for  r - 1 

(r-D(M-K) 

for  r > 2. 

4(2-p) 


(6.14) 


Note  that  from  relation  (3.5),  E[B2|F^=r]  can  be  obtained  using  the 
results  of  E[B|F^-r]  and  E[B2|F^-r]  from 

E[B|F^-r]  = E[B5^|F^=r]  + E[B2|F^=r]  for  all  r.  (6.15) 


6 . 4 Expected  Length  of  Times  related  to  the  Busy  Period 
Using  the  explicit  results  of  P[F_-r]  obtained  in  (4.18)  to  (4.19), 

a 

E[B|F^-r]  in  (6.8),  E[Bj^|F^-r]  in  (6.12)  and  E[B2|F^-r]  in  (6.14),  we 
can  obtain  the  expected  length  of  times  E[B] , E[C]  and  E[B£]  for  i-1,2 
related  to  the  busy  period. 

6.4.1  EfBI 

The  expected  length  of  the  busy  period,  E[B]  can  be  obtained  from 

00 

E[B]  - Z E[B[F^=r]  P[F^-r] . 
r-1 


(6.16) 
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Hence , we  have 

E[B]  - P[F^-l]E[Uo]  + (1  - P[F^=1]){E[U]^]  + E[U2]) 

+ {E[F^]  - l)E[V-^]  + {E[F^]  - 2 + P[F^=1])E[Y;lJ- 
or  equivalently 

E[B]  - P[D-0]E[Uol  + {1  - P[D-0]  )(E[U;l]  + E[U2]) 

+ E[D]E[V]^]  + {E[D]  - 1 + P[D-0]  )E[Yi^]  . (6.18) 


Substituting  P[F^  - 1]  and  E[F^]  in  (4.18)  to  (4.19)  and  (6.1)  into 
relation  (6.17),  or  P[D-1]  and  E[D]  in  (4.20)  to  (4.21)  and  (6.2)  into 
relation  (6.18)  yields 


E[B] 


E[Uq] 


(,«-N.,M) 

+ — {E[Uj^]  + E[U2l) 

(1-P^) 


+ 


(pM-N.,M) 

“o^ 


E[Vi] 


(1-P^)(1-P«-^) 


E[Y;^] 


for  p / 1 


(6.19) 

(M-N)  N N 

E[Uq]  + {E[Uj^l  + E[U2])  + E[V3_] 

M M M-K 


N K 

+ E[Y^]  for  p - 1. 

M(M-K) 


Finally,  substituting  E[Uq]  , E[U^^]  , E[U2],  E[Y^^]  and  E[V^]  in  (6.4)  to 
(6.8)  into  relation  (6.19)  and  simplifying,  we  obtain  the  same  E[B]  as 
derived  in  (2.34). 


6.4.2  E[C1 


From  relation  (2.30),  we  have 

E[C]  - N/A  + E[B] . 


(6.20) 


Thus,  substituting  E[B]  into  the  above  relation  yields  the  same  E[C] 


as 
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in  (2.35). 


6.4.3  E[B^-]  for  i-1 , 2 

and  £[62]  , the  expected  length  of  time  when  exactly  one 
service  station  is  operating  and  when  both  service  stations  are 
operating  simultaneously  during  the  busy  period  are  obtained  from 

00 

E[BJ  - E E[BjF^-r]  P[F^-r]  for  i-1,2.  (6.21) 

r=l 

Hence , we  have 

E[Bi]  - P[F^-1]  E[Uq]  + {1  - P[F^-l]HE[Uj^]  + E[U2]) 

+ (E[F^]  - 2 + P[F^-l])E[Yj^]  . (6.22) 

E[B2]  = { E[F^]  - 1 ).E[Vj^].  (6.23) 


Substituting  P[F^-1]  in  (4.18)  to  (4.19),  E[F^]  in  (6.1)  and  E[Uq] , 

E[Uj^]  , E[U2],  E[Yj^]  and  E[Vj^]  in  (6.4)  to  (6.8)  into  relations  in  (6.22) 
and  (6.23)  and  simplifying,  we  obtain  E[B^]  and  £[82]  as 


E[B.] 


1 

A 


Np(l-p^-^^)  - 

(l-p)(l-p^-‘^) 


for  p / 1 


N(M+K-N) 


(6.24) 


for  p - 1, 


E[B2] 


A(2-p)(l-/-*^) 


for  p / 1 


N 

A 


for  p “ 1. 


(6.25) 
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6 . 5 Steady- State  Probabilities 

Without  solving  all  of  the  steady- state  difference  equations  in 
(2.1)  to  (2.10),  it  is  possible  to  obtain: 

1)  the  steady- state  probability  that  i service  stations  are  operating 
during  the  busy  cycle,  denoted  by  P[0-i]  for  i-0,1,2, 

2)  the  steady-state  probability  that  no  service  stations  are  operating 
and  no  customers  are  in  the  system,  denoted  by  P(0,0), 

using  the  explicit  forms  of  E[I],  E[Bj^],  E[B2]  and  E[C]  . 

6.5.1  PfO-il  for  i=0.1.2 

The  probability  that  i service  stations  are  operating  during  the 
busy  cycle,  P[0-i]  can  be  expressed  by  the  ratio  between  the  expected 
length  of  time  when  i service  stations  are  operating  and  the  expected 
length  of  the  busy  cycle  such  as 

E[I] 

P[0-0]  = , (6.26) 

E[C] 

E[Bi] 

P[0-1]  - , (6.27) 

E[C] 

E[B2] 

P[0-2]  “ . (6.28) 

E[C] 

Thus,  substituting  E[I]  and  E[C]  into  (6.26)  and  simplifying,  we  obtain 
the  same  P[0-0]  as  in  (2.33).  However,  substituting  E[B^^]  , E[B2]  and 
E[C]  in  (6.24),  (6.25)  and  (2.35)  into  relation  (6.27)  and  (6.28)  yields 
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P[0-1]  - 


M+K-N 

for  p - 1, 

M+K-N+4 


for  p / 1 

(6.29) 


P[0-2] 


(l-p)p«-N+l(i-pN)(j^.K) 
N(2-p)(l-p^-l^)  - p^-N+l(l-pN)(^.K) 
2 

for  p - 1. 

M+K-N+4 


for  p 1 

(6.30) 


6.5.2  P(0,0) 


From  the  rate  diagram  of  the  system  appearing  in  Figure  2.1,  it  is 
easy  to  verify  that 


P[0  - 0]  - N P(0,0).  (6.31) 

Hence , 


P[0  - 0] 

P(0,0)  - . 

N 

Substituting  P[0-0]  in  (6.26)  into  relation  (6.32)  yields 

1 

P(0,0)  - . 

AE[C] 


(6.32) 


(6.33) 


Thus,  it  is  easy  to  verify  that  we  have  the  same  P(0,0)  as  in  (2.21)  if 
we  substitute  E[C]  into  relation  (6.33). 
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6 . 6 Expected  Number  of  Operating  Service  Stations 

Since  the  obtained  P[0=i]  (i=0,l,2)  is  the  probability  mass 
function  of  the  number  of  operating  service  stations,  the  expected 
number  of  operating  service  stations,  E[0]  is  obtained  from 

2 

E[0]  - 2 iP[0-i].  (6.34) 

i-0 

Substituting  P[0-i]  in  (6.26)  to  (6.28)  into  relation  (6.34),  we  obtain 

E[B^]  2 E[B2] 

E[0]  - + . (6.35) 

E[C]  E[C] 

Hence,  substituting  E[Bj^]  , E[B2]  and  E[C]  in  (6.24),  (6.25)  and  (2.35) 
into  relation  (6.35)  and  simplifying,  we  obtain  the  same  E[0]  as  in 
(2.25).  In  particular,  relation  (6.35)  will  be  used  to  construct  a 
total  expected  cost  function  in  Chapter  7. 

6 . 7 Expected  Number  of  Activations  or  Removals  of 

Service  stations  per  Unit  Time 

E[Ra^  ot  E[R^] , the  expected  number  of  activations  or  removals  of 
service  stations  per  unit  time  can  be  obtained  from  the  following  five 
methods  which  will  lead  the  same  E[R^^]  as  in  (2.38)  depending  on  the 
available  information. 

Method  1 : Using  E[F^]  and  E[C] , we  have 

E[R^]  - E[Rj.]  - . (6.36) 

E[C] 

Method  2 : Substituting  relation  (3.3)  into  (6.36)  yields 

E[D]  1 


E[RJ  = 


E[C] 


+ 


E[C] 


(6.37) 
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Note  that  relation  (6.37)  will  be  used  to  construct  a total  expected 
cost  function  per  unit  time  in  Chapter  7. 


Method  3 : From  relation  (3.3)  and  (6.23),  we  have 

E[B2] 

E[D]  - . (6.38) 

E[Vi] 

From  relation  (6.28),  we  also  have 

E[B2]  - E[C]  P[0  - 2] . (6.39) 

Therefore,  substituting  E[B2]  in  (6.23)  into  relation  (6.38)  yields 

E[C]  P[0  - 2] 

E[D]  - . (6.40) 

E[Vi^] 

Substitute  E[D]  in  (6.40)  into  relation  (6.37),  we  obtain 

P[0  - 2]  1 

E[R^]  - + . (6.41) 

^[^l]  E[C] 


Method  4 : 


Using  relation  (6.22),  (6.27)  and  (6.37),  we  have 

P[0-0]  1 P[F^-l]E[Uo] 

■ + . 

E[Yj^]  E[C]  E[C]E[Y3^] 

{1  - P[F^-1]){E[U3^]  ^f^2]  - E[Y^]] 

E[C]E[Yi] 


(6.42) 


Method  5 : From  relation  (6.18),  we  have  E[D]  as 

E[B]-E[U3^]+E[Yj^]-E[U2]+P[D-0]  {E[U3^]+E[U2]-E[Uq]-E[Y3^]  } 

E[D] (6.43) 

E[Y]^]+E[V^] 

Hence,  substituting  E[D]  in  (6.43)  into  relation  (6.37)  yields 

E[B]-E[U^]+2E[Y^]+E[V^]-E[U2]+P[D-0]{E[U^]+E[U2]-E[Uq]-E[Y]^]) 
E[R  ] 


E[C]{E[Y^]+E[V^]} 


(6.44) 
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In  the  next  chapter,  we  will  construct  a total  expected  cost 
function  using  some  of  the  system  characteristics  derived  in  this 
chapter  and  determine  the  optimal  customer  levels  (K*,N*,M*)  numerically 
based  on  minimizing  the  total  cost. 


CHAPTER  7 


ECONOMIC  ANALYSIS  OF  THE  M/M/2  SYSTEM 
OPERATING  UNDER  THE  TRIADIC  (0,K,N,M)  POLICY 

7 . 1 Introduction 

In  previous  chapters,  we  derived  several  important  system 
characteristics  for  predicting  system  behavior.  The  problem  that  we 
address  now  consists  in  determining  the  optimal  customer  levels 

^ ^ "fc 

(K  ,N  ,M  ) at  which  service  stations  should  be  adjusted  based  on  an 
economic  criterion.  The  common  criterion  for  the  determination  of  the 
optimal  operating  policy  in  the  area  of  multi -removable  queueing  models 
is  to  minimize  the  total  expected  cost  involved  in  system  operation. 
Several  types  of  cost  function  could  be  selected.  Usually,  the  total 
expected  cost  involves  as  components  the  expected  holding  cost  for 
waiting  customers,  the  expected  operating  cost  for  operating  service 
stations  and  the  expected  switching  cost  for  activations  and  removals  of 
service  stations. 

It  is  important  to  point  out  that  most  of  the  cost  functions  used 
in  the  area  of  multi -removable  queueing  models  are  constructed  under  two 
assumptions.  The  first  is  the  assumption  of  homogeneous  structure  for 
each  cost  element.  In  other  words,  the  same  amount  of  the  holding  cost 
for  each  customer  per  unit  time,  the  operating  cost  for  each  operating 
service  station  per  unit  time  and  the  switching  cost  for  each  activation 
or  removal  of  service  stations,  are  assumed  regardless  of  the  system 
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states.  The  other  involves  no  consideration  of  the  cost  for  service 
stations  removed  from  the  system  to  perform  auxiliary  tasks. 

In  this  chapter,  we  construct  a total  expected  cost  function  per 
unit  time  which  could  be  applicable  to  a variety  of  real  situations. 

The  cost  for  service  stations  removed  from  the  system  associated  with 
performance  of  auxiliary  task  is  included  in  the  cost  elements.  We 
assume  that  auxiliary  task  is  always  available  for  service  stations 
removed  from  the  system.  We  also  assume  that  the  structures  for  some  of 
cost  elements  vary  depending  on  the  system  states.  We  do  so  by  relaxing 
the  assumptions  of  homogeneous  structures  for  these  cost  elements. 

Following  the  construction  of  the  total  expected  cost  function  per 
unit  time,  we  determine  the  optimal  customer  levels  (K*,N*,M*) 
numerically  at  which  service  stations  should  be  adjusted  in  order  to 
minimize  the  total  expected  cost  function.  Finally,  we  perform  a 
sensitivity  analysis. 

7.2  Construction  of  the  Total  Expected  Cost  Function 

In  order  to  construct  the  total  expected  cost  function  per  unit 
time,  we  assume  that  the  cost  elements  consist  of  (1)  the  holding  cost 
for  waiting  customers,  (2)  the  operating  cost  for  operating  the  service 
stations,  (3)  the  activation  cost  and  the  removal  cost  for  activating 
and  removing  the  service  stations  and  (4)  the  cost  for  service  stations 
removed  from  the  system  associated  with  performance  of  auxiliary  tasks. 
We  also  assume  that  the  cost  for  service  stations  removed  from  the 
system  associated  with  performance  of  auxiliary  tasks  involves  (i)  the 
cost  to  start  performing  auxiliary  task  when  removed  from  the  system, 

(ii)  the  cost  to  perform  auxiliary  task  and  (iii)  the  cost  to  stop 
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performing  auxiliary  task  when  reactivated.  However,  the  structures  of 
all  cost  elements  except  the  holding  cost  vary  depending  on  the  number 
of  operating  service  stations. 

Next,  we  define  the  following  mathematical  representations  for  the 
cost  elements. 

(a)  a ; the  holding  cost  for  one  customer  in  the  system  per  unit  time  : 
The  holding  cost  can  be  considered  as  the  penalty  cost  for  delaying 
service  to  the  customers  who  are  waiting  in  the  system. 

(b)  (i— 1,2)  ; the  operating  cost  for  operating  the  service  stations 

per  unit  time  when  i service  stations  are  operating  : The  operating 

cost  is  incurred  by  the  operating  service  stations  to  provide  service 
for  the  customers. 

(c)  a ’’’ r ’ activation  cost  and  the  removal  cost  for  activating 

and  removing  the  service  stations,  respectively  : The  activation  cost 

is  incurred  each  time  when  a seirvice  station  starts  operating  while  both 
service  stations  remain  inoperative  temporarily.  The  removal  cost  is 
incurred  each  time  when  the  operating  service  station  is  removed  from 
the  system  while  the  other  service  station  is  already  removed  from  the 
system.  In  other  words,  the  activation  cost  is  incurred  once  at  the 
beginning  of  the  busy  period  and  the  removal  cost  is  incurred  once  at 
the  end  of  the  busy  period  during  the  entire  busy  cycle. 

(d)  ^ ; the  cost  to  start  performing  auxiliary  task  by  the  service 

stations  : This  cost  is  incurred  each  time  when  the  service  station 

removed  from  the  system  starts  performing  auxiliary  task  while  the  other 
service  station  is  still  operating.  In  other  words,  this  cost  is 
incurred  whenever  the  number  of  customers  in  the  system  reaches  K while 
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both  service  stations  are  operating  simultaneously.  Hence  this  cost  is 
incurred  D times  during  the  busy  cycle  where  D was  defined  in  Chapter  3. 

(e)  (i-1,2)  ; the  cost  to  perform  auxiliary  task  by  the  service 
stations  per  unit  time  when  i service  stations  are  performing  auxiliary 
task 

(f)  ^ ; the  cost  to  stop  performing  auxiliary  task  by  the  service 

stations  when  reactivated  : This  cost  is  incurred  each  time  when  the 

service  station  removed  from  the  system  stops  performing  auxiliary  task 
when  reactivated  while  the  other  service  station  is  already  operating. 

In  other  words,  this  cost  is  incurred  whenever  the  number  of  customers 
in  the  system  reaches  M while  one  service  station  is  operating.  Hence, 
this  cost  is  incurred  D times  during  the  busy  cycle. 

(g)  TC(K,N ,M)  ; the  total  expected  cost  function  per  unit  time  : This 

is  the  sum  of  (i)  the  expected  holding  cost  per  unit  time,  (ii)  the 
expected  operating  cost  per  unit  time,  (iii)  the  expected  activation 
cost  and  the  removal  cost  per  unit  time  and  (iv)  the  expected  cost  per 
unit  time  incurred  by  the  service  stations  removed  from  the  system 
associated  with  their  performance  of  auxiliary  tasks. 

Using  the  definitions  of  each  cost  element  above  and  the  system 
characteristics  obtained  in  the  previous  chapters,  the  total  expected 
cost  function  per  unit  time  can  be  represented  by 


TC(K,N,M)  - aE[J]  + 


1 


+ ^2 


E[C] 


E[C] 


E[C] 


E[I] 


E[D] 


(7.1) 


E[C] 


E[C] 


E[C] 
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where  E[J]  is  the  expected  number  of  customers  in  the  system, 

E[B^]  (i-1,2)  is  the  expected  length  of  time  when  i service 
stations  are  operating  during  the  busy  period, 

E[C]  is  the  expected  length  of  the  busy  cycle, 

E[I]  is  the  expected  length  of  the  idle  period, 

E[D]  is  the  expected  number  of  times  the  level  of  customers  in  the 

system  reaches  M while  one  service  station  is  operating 

during  the  busy  period. 


Let 


^ar  - ’’a  + ’’r- 


^st  - ^s  + 


(7.2) 

(7.3) 


Then,  substituting  E[J],  E[Bj^]  , E[B2],  E[I],  E[C]  and  E[D]  obtained  in 
(2.23)  to  (2.24),  (6.24),  (6.25),  (2.32),  (2.35)  and(6.2)  into  relation 
in  (7.1)  and  simplifying  yields  the  total  expected  cost  function  per 
unit  time  explicitly  as 


TC(K,N,M) 


ap(3-2p)  a(M+k+l) 

+ + 

(l-p)(2-p)  2 

aN(l-p^‘^){  (M+K-N)(2-p)  + 4 ) 

2{N(l-p^-^)(2-p)  - p^-N+l(l-pN)(M-K)) 

A(l-p)(2-p)[r^^(l-pM-K)  + 

N(l-p^-‘^)(2-p)  - p^-N+l(l-pN)(f^.K) 
(1-p)  (y32-;9i-5i)  (M-K)p^-N+l(l-pN) 

N(1-p^-K)(2-p)  - p^*N+l(l-/)(M-K) 
(1-p)(52-Si-/3i^)N(1-p^'^)(2-p) 
N(l-p^-^)(2-p)  - p^-N+l(l-pN)(M.K) 


+ 


for  p / 1, 
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M+3  6(N-1)  + (K+N)(K-N)  + M(N-4) 

2 3(M+K-N+4) 


TC(K,N,M)  - 


(;3^+fi^)(M+N-K)  + 2{^2^S2) 

+ 

(M+K-N+4) 

2A[r^j.(M-K)  + €g^N] 

+ for  p - 1. 

N(M-K) (M+K-N+4) 


(7.4) 


In  general , we  assume  that 

^2^  2^1-  (7.5) 

< S2  ^ 25]^.  (7.6) 

The  relations  above  are  similar  to  the  concepts  of  individual  and  Joint 
ordering  policies  in  inventory  theory  (see  Sivazlian  and  Stanfel[21] ) . 
However,  if  we  set 


- p2  - 2^, 

26-^  - 62  - 25, 
^ar  “ ^st  “ ^ • 


(7.7) 

(7.8) 

(7.9) 


then,  TC(K,N,M)  in  (7.1)  becomes 

TC(K,N,M)  - aE[J]  + /9E[0]  + 5{2-E[0])  + rE[R^] . (7.10) 

TC(K,N,M)  is  a common  form  for  the  total  expected  cost  function  used  in 
the  area  of  the  multi -removable  service  stations.  Note  that  each  cost 
element  maintains  the  homogeneous  structure.  However,  from  relation 
(2.25),  E[0]  is  not  a function  of  the  decision  variables  (K,N,M). 

Hence,  for  the  determination  of  optimal  customer  levels  which  minimize 
TC(K,N,M),  relation  (7.10)  is  equivalent  to 


TC(K,N,M)  - qE[J]  + rE[R  1 . 

a 


(7.11) 
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Next,  we  determine  the  optimal  customer  levels  (K*,N*,M*) 
numerically. 


7 . 3 Determination  of  the  Optimal  Customer  Levels 

As  shown  in  (7.4),  it  is  very  difficult  to  determine  the  optimal 

customer  levels  (K  ,N  ,M  ) analytically  which  minimize  the  total 

expected  cost  function  per  unit  time  because  of  the  complexity  of  the 

optimization  problem,  i.e.,  a nonlinear  integer  minimization  problem. 

^1* 

One  of  the  possible  way  to  obtain  (K  ,N  ,M  ) is  a numerical  approach. 

We  provide  three  numerical  examples  to  obtain  the  optimal  customer 
levels  for  given  specific  values  of  the  system  parameters,  A and  ^ and 
the  cost  elements,  a,  /9^(i-l,2),  5^(i-l,2),  and 

Example  1 : Determination  of  (K  ,N  ,M  ) given  the  value  of  p 

Suppose  p = 0.6.  Also  suppose  a=$1.0,  /9^=$3.0,  yS2“$5.0,  5]^-$1.0, 
52=$1.5,  t^j.=$10.0  and  eg^=$5.0. 

Table  7.1  exhibits  the  results  of  (K*,N*,M*)  and  its  corresponding 
TC(K*,N*,M*)  for  different  values  of  A and  p. 


Table 

7.1 

(K*,N*,M*)  given 

p = 0.6 

A 

4 

(K*,N*,M*) 

-ie  -ie  ’4e 

TC(K  ,N  ,M  ) 

3 

5 

(3.5,10) 

$ 8.7319 

6 

10 

(5,7,14) 

10.8573 

9 

15 

(7,8,18) 

12.4544 

15 

25 

(9,11,23) 

14.9422 

30 

50 

(13,15,32) 

19.4979 

Ill 


Example  2 : Determination  of  M*  given  K=N=1 

Suppose  A - 1.  Also  suppose  a=$1.0,  y0^-$3.O,  5^=$1.0, 

find  £^^=$5.0. 

Then,  M is  obtained  for  different  values  of  n and  is  shown  in 
Table  7.2. 


Table 

7.2  M* 

given  K-M-1 

A 

M 

M* 

TC(1,1,M*) 

1 

4 

13 

$ 9.0208 

1 

3 

11 

8.3889 

1 

2.5 

10 

7.9663 

1 

2.0 

8 

7.4883 

1 

1.5 

6 

7.0835 

1 

0.9 

4 

7.3003 

1 

0.7 

3 

7.9496 

1 

0.6 

3 

10.8437 

1 

0.5 

2 

15.6991 

Example  3 : Determination  of  (K*,N*,M*)  given  2S2-S2,  ^ar”^st 

Suppose  p - 0.6  where  A - 3.0  and  - 5.0.  Also  suppose  a-$1.0, 
y9j^-$3.0,  /92-$6.0,  63^-$!.  0,  52"$2.0  and  ’■ar"^st"^5  • ° • 

Then,  (K  ,N  ,M  ) and  TC(K*,N*,M*)  are  obtained  as  (1,1,9)  and 
$7.0573. 

As  shown  in  Table  7.1  and  7.2,  it  is  easy  to  see  that  the  optimal 
customer  levels  change  significantly  for  different  values  of  the  system 
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parameters,  A and  As  the  values  of  A and  n are  increased,  the 
optimal  customer  levels  are  also  increased. 

Next,  we  perform  a sensitivity  analysis. 

7.4  Sensitivity  Analysis 

Based  on  the  case  when  A - 3.0  and  - 5.0  in  Example  1,  we  perform 
a sensitivity  analysis  for  changes  in  the  optimal  customer  levels  along 
with  changes  in  one  of  the  following  cost  elements: 

(1)  P2’ 

(2)  52- 

(3)  -ar- 

(4)  €3,. 

The  values  of  the  other  cost  elements  remain  the  same  as  in  Example  1. 

The  values  of  (K*,N*,M*)  and  the  corresponding  TC(K*,N*,M*)  are 
shown  in  Table  7.3  to  7.6.  It  is  easy  to  see  that  the  optimal  customer 
levels  do  not  change  significantly  for  different  values  of  a specific 
cost  element.  However,  exact  variations  of  (K*,N*,M*)  are  not 
predictable . 

In  the  next  chapter,  we  will  discuss  topics  for  future  research  in 
the  area  of  the  controllable  M/M/2  queueing  system  operating  under  the 
triadic  (0,K,N,M)  policy. 
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Table  7.3  Sensitivity  analysis  with  changes  in  ^2 


^2 

(K*,N*,M*) 

TC(K*,N*,M*) 

$6.0 

(3,5,11) 

$8.7455 

5.5 

(3,5,11) 

8.7393 

5.0 

(3,5.10) 

8.7319 

4.5 

(3,5,10) 

8.7227 

4.0 

(3,5,10) 

8.7134 

3.5 

(2.5,10) 

8.7041 

3.0 

(2,5.10) 

8.6932 

Table  7.4 

Sensitivity  analysis 

with  changes 

52 

(K*,N*,M*) 

TC(K*,N*,M*) 

$1.0 

(3,5,10) 

$8.5227 

1.2 

(3,5,10) 

8.6064 

1.4 

(3.5,10) 

8.6901 

1.5 

(3,5,10) 

8.7319 

1.7 

(3.5,11) 

8.8156 

1.9 

(3.5,11) 

8.8982 

2.0 

(3,5,11) 

8.9393 
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Table  7.5 

Sensitivity  ana 

lysis  with  changes 

^ar 

"ie  'Ic  "ic 

(K  ,N  ,M  ) 

TC(K*,N*,M*) 

$10.0 

(3,5,10) 

$8.7319 

9.0 

(4,4,10) 

8.4439 

7.0 

(3,3,9) 

8.8059 

5.0 

(1,1,8) 

6.8417 

4.0 

(1,1,7) 

5.5942 

2.0 

(1,1,6) 

3.0329 

0.0 

(1,1,5) 

0.3737 

Table  7.6 

Sensitivity  analysis  with  changes 

^st 

•Je  'Ar  'Je 

(K*,N*,M*) 

TC(K*,N*,M*) 

$5.0 

(3,5,10) 

$8.7319 

4.0 

(3,5,10) 

8.7319 

3.0 

(3,5,9) 

8.7262 

2.0 

(3,5,8) 

8.7111 

1.0 

(3,5.7) 

8.6693 

0.0 


(4,5.6) 


8.5789 


CHAPTER  8 


EXTENSIONS  AND  TOPICS  FOR  FUTURE  RESEARCH 
8 . 1 Introduction 

Because  the  fluctuation  of  the  number  of  customers  in  the  system  in 
a Markovian  queueing  system  is  a random  walk,  the  properties  of  the 
gambler's  ruin  problem  can  provide  fundamental  information  to 
characterize  such  system.  Hence,  the  methodology  used  to  analyze  the 
busy  period  may  be  applicable  to  other  Markovian  queueing  models. 

In  this  chapter,  we  discuss  the  following  in  the  area  of  the 
controllable  M/M/2  queueing  system  operating  under  the  triadic  (0,K,N,M) 
policy  as  topics  for  future  research; 

(1)  development  of  efficient  algorithms  or  heuristic  methods  to  obtain 
the  optimal  customer  levels  which  minimize  the  total  expected  cost 
function  per  unit  time  derived  in  (7.4); 

(2)  analysis  of  the  model  with  a finite  system  capacity,  denoted  by  R; 

(3)  analysis  of  the  model  with  a finite  population  size,  denoted  by  Q. 
Besides  the  above,  construction  of  a total  cost  function  based  on  a 

discounted  cost  criterion,  analysis  of  the  model  with  priority  service 
discipline  and  analysis  of  the  model  with  system  dependent  service  rate 
are  also  suggested  as  topics  for  future  research. 


115 


116 


8 . 2 Development  of  Efficient  Algorithms  or  Heuristic  Methods 
to  obtain  the  Optimal  Customer  Levels 

As  discussed  in  Chapter  7,  it  is  very  difficult  to  determine  the 
optimal  customer  levels  (K*,N*,M*)  analytically  which  minimize  the  total 
expected  cost  function  TC(K,N,M)  because  of  the  complexity  of  the 
problem.  Hence,  one  of  the  possible  ways  to  obtain  (K*,N*,M*)  is  a 
numerical  approach.  However,  the  numerical  approach  necessitates  a huge 
amount  of  computational  effort  because  the  values  of  TC(K,N,M)  for  all 
possible  values  of  the  decision  variables  (K,N,M)  have  to  be  evaluated. 
This  situation  leads  towards  the  development  of  efficient  algorithms  or 
heuristic  methods  to  obtain  the  optimal  customer  levels. 

From  relation  (7.1),  TC(K,N,M)  is  given  by 

E[B^]  E[B2]  1 

TC(K,N,M)  = oE[J]  + + ^2 + ’’ar 

E[C]  E[C]  E[C] 

E[Bj_]  E[I]  E[D] 

+ 6]^ + 5 2 + . (8.1) 

E[C]  E[C]  E[C] 

It  is  clear  that  TC(K,N,M)  depends  of  the  values  of  the  system 
parameters  and  the  cost  elements.  Hence,  an  analysis  of  the  functional 
behavior  of  E[J],  E[B^]/E[C],  E[B2]/E[C] , 1/E[C] , E[I]/E[C]  and 
E[D]/E[C]  might  provide  fundamental  information  to  predict  the  behavior 
of  TC(K,N,M).  Also  it  might  be  necessary  to  analyze  which  cost  elements 
are  critical  in  evaluating  TC(K,N,M)  for  different  values  of  the  system 
parameters.  It  may  be  possible  to  develop  efficient  algorithmic 
procedures  or  heuristic  methods  to  obtain  the  optimal  customer  levels  if 
additional  analysis  is  performed  in  the  structure  of  the  cost  function 
which  takes  the  form  of  a nonlinear  integer  optimization  problem. 
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8 . 3 Analysis  of  the  Model  with  a Finite  System  Capacity 

Suppose  that  the  system  capacity  is  limited  to  R which  is  a 
positive  integer  and  is  greater  than  M.  Then,  no  customers  are  allowed 
to  enter  the  system  when  R customers  are  in  the  system.  As  shown  in  the 
rate  diagram  of  the  model  appearing  in  Figure  8.1,  customers'  mean 
arrival  rate  is  still  independent  of  the  system  states  because  of  an 
infinite  population  size.  The  only  difference  between  this  model  and 
the  model  with  an  infinite  system  capacity  and  an  infinite  population 
size  analyzed  in  the  previous  chapters  is  the  fluctuation  of  the  number 
of  customers  in  the  system  when  both  service  stations  are  operating 
simultaneously.  This  implies  that  the  results  related  to  the 
fluctuation  of  the  number  of  customers  in  the  system  when  one  service 
station  is  operating  obtained  in  the  previous  chapters  are  directly 
applicable  to  this  model.  Hence,  this  model  can  be  easily  analyzed  by 
deriving  necessary  information  for  the  fluctuation  of  the  number  of 
customers  in  the  system  when  both  service  stations  are  operating 
simultaneously. 

The  steady- state  behavior  of  the  model  can  be  studied  by  solving 
the  difference  equations  for  the  steady- state  joint  probability  mass 
function  P(i,j)  for  i-0,1,2  and  j-0 , 1 , 2 , • • • , R.  The  difference  equations 
consist  of  (2.1)  to  (2.8),  (2.9)  and 

(A+2m)P(2,J)  - 2MP(2,j+l)  + AP(2,j-l)  (K+2<j^-l,  R^M)  , (8.2) 

2mP(2,R)  = AP(2,R-1) . (8.3) 

The  values  of  p and  q which  are  the  elements  of  the  one  step 
transition  probability  matrix  of  the  number  of  customers  in  the  system 
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number  of  operating  service  stations 


2 10 


Figure  8.1  The  rate  diagram  of  the  model 
with  a finite  system  capacity 
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when  both  service  stations  are  operating  simultaneously  remain  the  same 
as  in  (3.11).  However,  the  one  step  transition  probability  that  the 
number  of  customers  in  the  system  will  reach  R-1  after  it  reaches  R is 
given  by 

P[J(n+l)-R-l|j(n)-R]  - 1 for  all  n.  (8.4) 

Hence,  the  one  step  transition  probability  matrix  of  the  number  of  the 
customers  in  the  system  when  both  service  stations  are  operating 
simultaneously  has  a reflecting  barrier  at  the  state  R as  shown  in 
Figure  8.2. 

Under  steady-state  conditions,  P(Vj^)  which  is  the  probability  that 
the  number  of  customers  in  the  system  will  decreased  to  K after  it 
reaches  M when  one  service  station  is  operating,  is  one.  Therefore,  the 
probability  mass  function  of  the  number  of  activations  or  removals  of 
service  stations  during  the  busy  period  is  exactly  the  same  as  P[F^— r] 
obtained  in  (4.18)  and  (4.19). 

In  order  to  obtain  the  Laplace  transform  of  the  probability  density 
function  of  the  busy  period,  the  probability  generating  function  of  the 
number  of  transitions  during  G2(M,z)  has  to  be  recomputed.  This 
will  lead  to  deriving  the  Laplace  transform  of  the  probability  density 
function  of  62(111, z)  also  satisfies  a second  order  difference 

equation  with  boundary  condition  given  in  (5.29)  and  (5.30).  Thus, 

G2(M,z)  “ pz  G2(M+1,z)  + qz  G2(M-l,z)  for  K<M<R-1  (8.5) 

with  boundary  condition 

G2(M,z)  - 1.  (8.6) 

However,  another  boundary  condition  related  to  G2(R,z)  has  to  be 
evaluated  based  on  the  one  step  transition  probability  matrix  of  the 


K 

1 

q 

0 

0 

0 

0 

0 

0 

0 
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K+1  K+2 
0 0 


0 

q 


p 

0 


0 0 
0 0 
0 0 

0 0 
0 0 
0 0 


M-1  M M+1 

0 0 0 

0 0 0 

0 0 0 


q 

0 


p 

0 


p 

0 


0 0 0 
0 0 0 
0 0 0 


R-2  R-1  R 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 

0 0 0 


p 

q 

0 


0 0 

0 p 

1 0 


The  one  step  transition  probability  matrix  of 
the  number  of  customers  in  the  system  when  both 
service  stations  are  operating  where  the  values 
of  p and  q are  shown  in  (3.11) 
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number  of  customers  in  the  system  in  which  the  state  R forms  a 
reflecting  barrier.  Thus,  the  difference  equation  (8.5)  is  solvable 
directly.  Substituting  (>i(2+p) }/{ s+/i(2+p) ) instead  for  z in  the 
solution  G2(M,z)  yields  the  Laplace  transform  of  the  probability  density 
function  of  Therefore,  substituting  the  Laplace  transform  of  the 

probability  density  function  of  and  the  existing  Laplace  transforms 
of  the  probability  density  functions  of  Uq,  U2  and  in  (5.65)  to 
(5.69)  into  relation  (5.64)  yields  the  Laplace  transform  of  the 
probability  density  function  of  the  busy  period.  Note  that  setting  K-1, 
N-1  and  M-2  yields  the  Laplace  transform  of  the  probability  density 
function  of  the  busy  period  for  the  ordinary  M/M/2  queueing  system  with 
a finite  system  capacity  R. 

The  total  expected  cost  function  per  unit  time  can  also  be  obtained 
by  recomputing  E[J],  £[82]  and  E[C] , and  then  substituting  them  into 
relation  (7.1). 

8 . 4 Analysis  of  the  Model  with  a Finite  Population  Size 

Suppose  that  the  population  size  or  the  total  number  of  possible 
customers  is  Q which  is  a positive  integer  and  is  greater  than  M.  As 
shown  in  the  rate  diagram  of  the  model  appearing  in  Figure  8.3, 
customers'  mean  arrival  rate  depends  on  the  number  of  customers  in  the 
system. 

For  the  steady- state  behavior  of  the  model,  the  joint  probability 
mass  function  P(i,j)  for  i-0 , 1 , 2 and  j-0 , 1 , 2 , • • • ,Q  has  to  be  recomputed 
using  the  following  difference  equations: 
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Figure  8.3  The  rate  diagram  of  the  model 
with  a finite  population  size 
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(Q-j)AP(i,j)  =.  (Q-j-l)AP(i,j+l)  (i=0,  0<j^-2),  (8.7) 

QAP(i,j)  - (i-0,  j-0),  (8.8) 

(M+(Q-j)A)P(i,j)  - MP(i,j+l)  (i-1,  j=l),  (8.9) 

{/i+(Q-j)A)P(i.j)  - (Q-j+l)AP(i,j-l) 

+ /iP(i,j+l)  (i-1,  2<j^-2,  j/K,M),  (8.10) 

{/i+(Q-j)A)P(i,j)  - (Q-j+l)AP(i,j-l)  + MP(i,j+l) 

+ 2/iP(i+l,j+l)  (i-l.j-K),  (8.11) 

{M+(Q-j)A)P(i,j)  “ (Q-j+l)AP(i,j-l)  + MP(i,j+l) 

+ (Q-j+l)AP(i-l,j-l)  (i-l,j-N),  (8.12) 

{/i+(Q-j)A)P(i,j)  - (Q-j+l)AP(i,j-l)  (i-l.j-M-1),  (8.13) 

{2M+(Q-j)A)P(i,j)  - 2MP(i,j-l)  (i-2,j-K+l),  (8.14) 

{2M+(Q-j)A}P(i,j)  - (Q-j+l)AP(i,j-l)  + 2/iP(i,j-l) 

(i=2,  K+2<j<Q-l,  j^M),  (8.15) 

{2/i+(Q-j)A}P(i,j)  - (Q-j+l)AP(i,j-l)  + (Q-j+l)AP(i-l,j-l) 

+ 2/iP(i,j+l)  (i-2,j-M),  (8.16) 

2MP(i,j)  - AP(i,j-l)  (i=2,j=Q).  (8.17) 


In  order  to  analyze  the  busy  period,  suppose  that  i (i-1, 2)  service 
stations  are  operating.  Let  for  j=l,2,---,Q  be  the  random  variable 
which  represents  the  length  of  time  from  the  instant  when  the  number  of 
customers  in  the  system  reaches  j to  the  instant  when  the  next 
transition  occurs.  Then,  the  T^^j's  are  independently  and  exponentially 
distributed.  It  can  be  verified  that  the  parameters  for  the  probability 
density  functions  of  the  ' s are  (Q-j)A  + i/i.  Let 

Pj^^  - P[J(n+l)-k|  J(n)=j  ] for  all  n.  (8.18) 

In  order  to  characterize  the  fluctuation  of  the  number  of  customers  in 


the  system,  the  following  one  step  transition  probabilities  of  the 
number  of  the  customers  in  the  system  can  be  written  down: 
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for  0=1  and  j=l , 2 , • • ,Q-1 , (8.19) 

for  0-2  and  j-1 , 2 , • • ,Q-1 , (8.20) 

for  0—1  and  j=l,2, • • ,Q,  (8.21) 

for  0-2  and  j-l,2,--,Q.  (8.22) 

The  fluctuation  of  the  number  of  customers  in  the  system  during  Uq, 

U2  and  are  also  dictated  by  the  one  step  transition  probability 
matrix  appearing  in  Figure  3.4.  However,  the  values  of  p and  q have  to 
be  replaced  by  pj  and  pj  respectively,  in  (8.19)  and  (8.21). 

The  one  step  transition  matrix  of  the  number  of  customers  in  the  system 
during  can  also  be  obtained  using  the  one  step  transition  probability 
matrix  shown  in  Figure  8.1  by  replacing  p and  q by  pj  and  pj  in 

(8.20)  and  (8.22)  after  changing  R to  Q. 

The  probability  mass  function  of  the  number  of  activations  or 
removals  of  service  stations  during  the  busy  period  can  be  obtained  from 
relation  (4.7)  because  of 

P(Vj^)  - 1.  (8.23) 

However,  the  values  of  P(Uq)  , P(U^)  , P(U2)  and  P(Yj^)  have  to  be 
recomputed.  For  example,  P(Uq)  is  obtained  by  solving  the  following 
second  order  difference  equation  where  'I>q(N)  is  equal  to  P(Uq)  : 


(Q-j)A 
(Q-j)A  + M 
(Q-j)A 
(Q-j)A  + 2m 


Pj.j-1 


(Q-j)A  + M 
2m 

(Q-j)A  + 2m 
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^>0<N)  - Pn,N+1  Pn,N-1  1<N:^-1 


(8.24) 


with  boundary  conditions 


$o(0)  = 1 

$q(M)  - 0. 


(8.25) 


(8.26) 


The  values  of  pjj  Pn  N-1  shown  in  (8.19)  and  (8.21). 

The  Laplace  transform  of  the  probability  density  function  of  the 
busy  period  is  obtained  from  relation  (5.63)  as 


Then,  setting  K-1,  N-1  and  M-2  in  fg(s)  yields  the  Laplace  transform  of 
the  probability  density  function  of  the  busy  period  for  the  ordinary 
M/M/2  queueing  system  with  a finite  population  size  Q.  However,  a new 
methodology  has  to  be  developed  in  order  to  derive  the  explicit  forms  of 
the  Laplace  transform  of  the  probability  density  functions  of  Uq, 

^2’  ^1  ^1'  Since  the  T^j's  are  not  identically  distributed, 

the  methodology  used  in  Chapter  5 to  obtain 
fy^(s)  and  fy^(s)  is  not  applicable  to  this  model. 
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APPENDIX  A 


DERIVATION  OF  THE  DISTRIBUTION  FUNCTION  OF  B,  B^^  AND  B2 


In  this  appendix,  we  provide  some  analytic  results  for  the 
distribution  function  of  B,  B^,  and  B2  discussed  in  Subsection  5.4.2. 

Inverting  the  probability  generating  functions  of  Lq,  L^,  L2 , 
and  obtained  in  (5.19)  to  (5.22)  and  (5.35)  yields  the  probability 
mass  functions  of  the  number  of  transitions  during  Uq,  U^^,  U2 , and 
for  m=l,2,---  as  (see  e.g.  Feller[10]) 

(l-p^)2"'p("'-N)/2  M-1 

P[LQ=m] — — E sin(j7r/M)sin(jNVM)cos®'^(j7r/M)  , (A.l) 

(1-P^-N)M(1+p)“>  j-1 


P[L2=ra] 


(1-P^)2“p^'"+^-N)/2  m-1 

— — — S sin(jVM)sin{j  (M-N)7r/M}cos'”'^(j7r/M)  , 

(p^-N.pM)M(l+p)“  j=l 


(A. 2) 

(1-P^)2^("'-‘^)/2  m-1 

— — S sin(jVM)sin(jKVM)cos"''^(j7r/M)  , (A. 3) 

(l-p"-‘^)M(l-hp)“  j-1 


P[L3-m] 


P[L^-m] 


(l-p^)2®p(‘"+^-^)/2  M-1 


M-K  M 


(p----p-)M(l+p)“  j-1 


2 sin(j7r/M)sin{  j (M-K)7r/M)cos“*'^(j7r/M)  , 


M-K 


m 


/ 


m 


\ (m-M+K)/22(m-t-M-K)/2 


(m+M-K)/2 


(A. 4) 
(A. 5) 


(2+p)“ 

The  length  of  times  between  two  successive  transitions  during  the 
busy  period  are  identically,  independently  and  exponentially  distributed 
with  parameter  (X+fi)  when  one  service  station  is  operating  and  with 
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parameter  (A+2/i)  when  both  service  stations  are  operating 
simultaneously.  Hence,  fy^(x|LQ=m)  has  a gamma  distribution  with 
parameter  (A+^)  and  of  order  m such  that 


(x+n) 

r(m) 


m 


j^m-lg-(A+M)x 


for  0<x<®. 


(A. 6) 


Substituting  P[LQ-m]  and  fy^(x|LQ=m)  derived  in  (A.l)  and  (A. 6)  into 
relation  (5.49)  which  is  given  by 


^ P[Lo-™]. 

m— 0 


(A. 7) 


we  obtain  fy^(x)  for  0<x<«>  as 


= 


e-(A+/^)X(i./)  « {2(A+m)x)“/“-N)/2  m-1 

^ 2:  0o(j.m) 

(l-p”''^)Mx  m=l  r(m)(l+p)“  j=l 


(A. 8) 


where 


-1/. 


“ sin(j5T/M)sin(jNVM)cos  ’ (j^r/M) 


Similarly,  we  obtain 


fu^(x) 


^-(A+/i)X(i.^M)  „ {2(A+m)x)'"/"'+^-N)/2  m-1 

S S 01  (j  ,m) 

r(m)(l+p)'"  j=l 


(A. 9) 


where 


m-1/ .. 


0]^(j,m)  - sin(j7T/M)sin(  j (M-N)7t/M)cos  ' (jVM). 


^-(A+M)X(i.pM)  „ {2(A+M)x}'"p('”-K)/2  M-1 

M Vf  ^ ni  ^ ®2^^ 

(l-p“'^)Mx  m-1  r(m)(l+p)®  j-1 


(A.  10) 


where 


m-1/4. 


e2(j,m)  - sin(j7T/M)sin(jK7r/M)cos  ’ (jVM), 
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fy^(x) 


S S G^Cj.m) 


(p^'^-p^)Mx  m=l  r(m)(l+p)‘"  j=l 


(A. 11) 


where 

03(j,m)  - sin(j7r/M)sin{j  (M-K)7r/M)cos"*'^(j?r/M)  , 

« / m \ /n‘-M+K)/22(m+M-K)/2(;^^2M)V-l 

f„  (x)  - (M-K)e-^^+2M)x  2 . 

1 ra-0  \(m+M-K)/2/  mr(m) (2+^)“ 

^ ' (A. 12) 

Using  relations  in  (5.59)  to  (5.62)  yield  the  probability  density 

function  of  B,  fg(x)  for  0<x<®  as 

fg(x)  - fuQ(x)P[F^-l]  + 

00 

S fu^(x)*f^^*(’"-^)(x)*fY^*(’'-2>(x)*fu2(x)  P[F^-r].  (A. 13) 

From  relation  (3.6),  the  probability  density  functions  of  for 
0<x<®  is  represented  by 

00 

fg^(x)  - fuQ(x)P[F^-l]  + fy^(x)*fYj^*^’'-2)(x)*fu2(x)  P[F^-r]  , (A. 14) 

Since  B2  has  a concentration  of  mass  at  x-0,  we  obtain  from 
relation  (3.7) 


P[F^-1] 


for  x-0 


2 f *(r-l)(x)  P[F^-r] 
r-2  ^ 


(A. 15) 


for  0<x<®. 
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